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ABSTRACT 

We present the first results of a new Keck spectroscopic survey of UV faint Ly- 
man break galaxies in the redshift range 3 < z < 7. Combined with earlier Keck and 
published ESO VLT data, our spectroscopic sample contains more than 600 dropouts 
offering new insight into the nature of sub-L* sources typical of those likely to dom- 
inate the cosmic reionisation process. In this first paper in a series discussing these 
observations, we characterise the fraction of strong Lya emitters within the continuum- 
selected dropout population. By quantifying how the "Lya fraction" varies with red- 
shift, we seek to constrain changes in Lya transmission associated with reionisation. In 
order to distinguish the effects of reionisation from other factors which affect the Lya 
fraction (e.g. dust, ISM kinematics), we study the luminosity and redshift-dependence 
of the Lya fraction over 3 < z < 6, when the IGM is known to be ionised. These results 
reveal that low luminosity galaxies show strong Lya emission much more frequently 
than luminous systems, and that at fixed luminosity, the prevalence of strong Lyman-a 
emission increases moderately with redshift over 3 < z < 6. Based on the striking cor- 
relation between blue UV slopes and strong Lya emitting galaxies in our dataset, we 
argue that the Lya fraction trends are governed by redshift and luminosity-dependent 
variations in the dust obscuration, with likely additional contributions from trends in 
the kinematics and covering fraction of neutral hydrogen. We find a tentative decrease 
in the Lya fraction at z ~ 7 based on the limited IR spectroscopic data for candidate 
z ~ 7 galaxies, a result which, if confirmed with future surveys, would suggest an 
increase in the neutral fraction by this epoch. Given the abundant supply of z and 
Y-drops now available from deep Hubble WFC3/IR surveys, we show it will soon be 
possible to significantly improve estimates of the Lya fraction using optical and near- 
infrared multi-object spectrographs, thereby extending the study conducted in this 
paper to 7 < z < 8. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: high-redshift 



1 INTRODUCTION 

Considerable observational progress has been achieved in re- 
cent years in the study of star-forming galaxies seen beyond 
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z ~ 3, a period corresponding to ~ 2 Gyr after the Big 
Bang. It is now clear that this is a period of rapid galaxy 
evolution and a number of key results have emerged from 
recent multi-wavelength surveys. 

For the colour-selected z > 3 Lyman break galaxies 
(LBGs), it is now established from various independent sur- 
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veys that the star formation density, deduced from rest- 
frame UV luminosities, declines monotonically w ith red- 
shift (e.g. IStanwav et~atl 120031 : iBunker et all 120041 ) largely 
as a result of a corr esponding fading of the characteristic 
UV luminosity (e.g. [Ouchi et al.ll2004l;lYoshida et al.ll2006l ; 



iBouwens et~atll2006i . l2007l : iMcLure et al.ll201(F The associ- 
ated stellar mass density in Lyman break galaxies, deduced 
from near-infrar ed Spitzer photometry, increases by ~1 dex 
from z ~ 6 to 4 (|Evles et al.ll2007l ; IStark et alj|2009h . As the 
rate of change of stellar mass is governed by ongoing star 
formation, it is useful to relate the two measures and such a 
comparison indicates a rapid duty-cycle of star formation ac- 
tivity at this time, unlike the mor e continuous mod es seen for 
equivalent sources at z ~ 2 — 3 l|Stark et alj|2009h . By con- 
trast, the redshift-dependent luminosity function of narrow- 
band selected Lyman alpha emitters, s hows no equivalen t 
decline with redshift over 3 < z < 6 l|Ouchi et al.| [2008). 
suggesting an increasing fraction of line emitters amongst 
the star forming population at early times. Moreover, de- 
tailed studies of the slope of the UV continuum in z > 3 
Lyman break gal axies indicates a decreasing dust content at 
earlier times (e.g. IStanwav et al.ll2005l ; IBouwens et al.ll200q . 
2009bl) as well as a lu minosity dependence at z ~ 3 (e.j 



Reddv fc St cidcl 200^). Conceivably the combination of a 
reduced dust content and a shift to more intense, shorter- 
term star formation at high redshift, can explain these var- 
ious redshift-dependent trends. 

Notwithstanding this considerable progress, a major 
concern is that the above conclusions rest largely on de- 
ductions made with photometric data, particularly for the 
Lyman break population. Quite apart from the possibility 
of low redshift interlopers lying within the photometric sam- 
ples (a problem that increases for drop-out selected sam- 
ples at redder wavelengths), the wanted physical measures 
of the star formation rate, dust content and stellar mass 
are all rendered uncertain by the absence of precise spectro- 
scopic redshifts. While considerable effort has been invested 
in th e spectroscopic study of z ~ 3 Lyman b reak galaxies 
(e.g. IShaplev et all 120031 ; ISteidel et all 120031 ; iQuider et ail 
l20ld ). comparably little spectroscopy has b een achieved for 
higher redshift samples. ISteidel et al.l |T999) obtained spec- 
troscopic redshifts for nearly 50 bright (I < 25) z ~ 4 
LBGs. Most surveys of z ~ 5 — 6 V and i'-drops have gen- 
erally involved relatively small samples, typically compris- 
ing fewer than 10 sources (e.g. IStanwav et al.l |2004| . 120071 ; 
Undo et al.ll2007l ; iDow-Hygelund et al.l 120071 ) Recent deep 
HST ACS Grism observations of the Hubble Ultra Deep 
Field have all owed the spectra of faint z — 5 LBGs to be 
characterised (|Rhoads et al J 12009). albeit at very low spec- 
tral resolution, resulting in 39 redshift confirmations. The 
largest spectroscopic sample of 4 < z < 6 LBGs thus far 
published is the VLT/FORS2 survey of the Chandra Dee p 
Field South (|Vanzella et al.ll2002l . |2005| . l200d . 120081 . 120091) . 
This survey represents a major step forward, targetting 195 
B, V and i'-drop galaxies and securing high redshifts for 99 
of them. In addition, a recent campaign with the VIMOS 
spectrograph on the VLT has also targeted bright z > 3.5 
LBGs (i 775 < 25) in the Chandra Deep Field South, con- 
firming redshifts for 20 bright sources at 3.5 < z < 5.5 
jBalestra et al.ll20ld ). 

We build on these strides in this paper, the first in a se- 
ries presenting the results of a new Keck survey of 3 < z < 7 



LBGs selected photometrically in the GOODS fields. The 
overall goal is to improve our understanding of the evolution 
of star-forming galaxies during the first 2 Gyr of cosmic his- 
tory. Fully exploiting the 10 metre Keck aperture, we have 
designed our survey to target intrinsically fainter sources 
than those reached in the VLT/FORS2 survey, thereby com- 
plementing that effort. Spectroscopy spanning a wide range 
of intrinsic luminosities is very important if we seek to under- 
stand earlier e xamples of the lumino sity-dependent trends 
seen at z ~ 3 (|Reddv fc Steidell 120091 ). As discussed below, 
it is equally important to target and study sub-luminous 
star forming sources at early times, as these may be typ- 
ical of those galaxies responsible for cosmic reionisation 
dBouwens et al.ll2007l; lOuchi et al.ll2009l ; IBunker et al.ll20ld ; 
lOesch et al.ll20ld ; iMcLure et al.ll20ldT 

Understanding evolution in the demographic trends of 
star-forming galaxies over 3 < z < 7 is vital if we are to use 
galaxies as tracers of cosmic reionisation, a cosmic event 
which is now a frontier of observational cosmology. Knowl- 
edge of when reionisation occurred is crucial to our under- 
standing of the nature of the earliest UV-emitting sources 
as well as the discrepancy between the observed number of 
dwarf galaxie s and those expected from cosmological sim- 
ulations (e.g. ISalvadori fc Ferraral 120091 ). Current observa- 
tional constraints are limited. WMAP measurements indi- 
cate that the uni verse may have been partially ionised a s 
early as z ~ 11 jDunklev et al l 120091 : lLarson et ai]|201Ch . 
while observations of transmission in the spectra of quasars 
reveal that intergalactic hydrogen must be highly ionised 
below z < 6. While the discovery of an accelerated decline 
and increased variance in the mean transmitted flux from 
quasars initially led many to suggest that the in tergalactic 
medi um (IGM) is still partially neutral at z ~ 6.2 (|Fan et al.l 
2006), recent work has demonstated that these results do 
not necessarily re quire a sudden cha nge in the ionisation 
state of the IGM (|Becker et al.ll2007l ). Given the difficulty 
in locating quasars at z > 7, it seems unlikely that quasar 
spectroscopy will constrain the epoch when the bulk of the 
IGM was reionised in the near future. 

Lya emitting galaxies o ffer a valuable add i tiona l 

probe of the IG M [e .g.. iRhoads fc Malhotral 1200 ll : 

iMalhotra fc Rhoadsl [2004 : iKashikawa et al.l l2006h ~ In prin- 



ciple, the test is straightforward to apply. Young galax- 
ies emit copious amounts of Lya photons, which are 
resonantly scattered by neutral hydrogen. Hence as we 
probe the regime when the IGM becomes significantly neu- 
tral, the fraction of star-forming galaxies showing strong 
Lya emission should decrease ( e.g. Haiman fc Spaans 1999; 
I Santosl |2004 iFurlanetto et all l200d: iMcQuinn et alj |2007l: 
i Mesinger fc Furlanettd |200S| ; llliev et all 120081 : iDaval et all 
Recent measurements of the luminosity function of 
Lya emitters (LAEs) selected via narrowband imaging have 
reve aled a tantalising decline between z = 5.7 and z = 7. 
(e.g. IKashikawa et alj|200d : live et alj|2005 [ Ota et al ] |2008l ). 
offering possible evidence that the ionisation state of the 
IGM evolves over 6 < z < 7. 

Regardless of the validity of claims for an increasing 
fraction of neutral hydrogen over 6 < z < 7, in prac- 
tice the interpretation of the Lya test is more complex. 
While the ionisation state of the IGM affects the Lya LF, 
so does evolution of a multitude of oth er properties in- 
trinsic to the sampled population (e.g., IVerhamme et all 
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l200rj, 120081; lAtek et alj|2008h . Evolution in the dust content 
(Bouwcns ct al . 2009b), the column density, kinematics, and 
geometrical distribution (generally de scribed as the "cov- 
ering fraction") of neut r al hydrogen ( Shaplev et al.l 120031 ; 
lOuider et al]|200Sl l20ld ; ISteidel et aljboioh can "each play 



a key role. Also important for the transmission of Lya pho- 
tons is the density of the IGM, which evolves continuously 
with redshift, and the stellar initial mass function of galax- 
ies, for which few robust constraints exist at high-redshift. 

The existence of these complicating factors highlights 
the importance of understanding how the prevalence of Lya- 
emitting galaxies varies just after reionisation in addition to 
characterising the decline that may occur during reionisa- 
tion itself. With this goal in mind, we seek to construct an 
independent measure of the redshift evolution of Lya emit- 
ting galaxies over a large redshift b aseline, complementary 
to past efforts (e.g. lOta et al"1l2008l ). Here we introduce the 
method of measuring the fraction of strong Lya emitters 
(hereafter the "Lya fraction") within the LBG population. 
Applying the Lya fraction test to a large sample of LBGs 
has many advantages. Firstly, the LBG samples are already 
in p lace over 4 < z < 8 owing to deep surveys with HST 
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2010b); hence the only time investment required is follow-up 
spectroscopy. This spectroscopy not only provides a sample 
of Lya emitters with known spectroscopic redshifts, but it 
also provides information on the kinematics of the ISM (e.g . 
IShaplev et aT1l2003l ; IVanzella et ai1l2009l ; ISteidel et al]|201Ch 
for individual bright sources (and for faint systems via com- 
posite spectra) and improves estimates of dust obscuration 
via UV colors, both of which are key factors governing the 
transmission of Lya photons. By improving our understand- 
ing of how these properties change with time, we can begin 
to isolate the effect of reionization on the evolution in the 
Lya fraction. Additionally, the Lya fraction is insensitive 
to the declining number density of star-forming galaxies, 
in contrast to the Lya luminosity function test (which re- 
quires comparison to the UV luminosity function to account 
for this degeneracy). Naturally, the Lya fraction test has 
its own complications, but we show that these can be cor- 
rected for (see §3.3 and 3.4), and we thus argue that this 
test will provide valuable constraints on reionization as new 
7 < 2 < 10 LBGs emerge in the next several years. 

Our goal is therefore to obtain a robust measure of the 
luminosity and redshift dependence of the Lya fraction when 
the IGM is highly-ionised (z ~ 3 — 6) and to compare it to 
the Lya fraction at progressively earlier times. If the IGM 
ionisation state evolves significantly at z > 6 as possibly im- 
plied by the narrowband Lya results, we would expect the 
measured fraction of LBGs with strong Lya emission to be 
lower at z ~ 6 — 8 than expected from extrapolating the 
trends seen over 3 < z < 6. By placing the evolution of the 
Lya fraction at 3 < z < 6 in the context of the evolution 
of the well-characterised LBG parent population (e.g. vari- 
ation of dust extinction with redshift and luminosity), we 
will calibrate the relative importance of factors other than 
reionisation on the transmission of Lya photons. 

Our Lya fraction test is ideally suited to our large spec- 
troscopic sample of continuum-selected LBGs spanning the 
redshift range 3 < z < 7. The current Keck sample (includ- 
ing a sample of archival i'-drops to be presented in Bunker et 
al. 2010, in preparation) consists of 455 B, V, i' , and z-drops 



(photometrically-selected to lie at 3.5 < z < 7.0) spanning 
a wide range in UV luminosity (to Muv — —18). We com- 
bine this sample with a more luminous publically-available 
FORS2 dataset of 195 sources satisfying photometric crite- 
ria similar to those adopted for the Keck surveys. Using this 
large spectroscopic database, we address the primary goal of 
this first paper in our series - to compute the Lya fraction 
as a function of luminosity and redshift. This will enable us 
to identify the principal factors governing its evolution over 
3 < 2 < 6, prior to extending the test to 6 < z < 8, where 
we can probe changes in the IGM. 

The plan of the paper is as follows. In §2, we describe the 
target selection and spectroscopic observations, and present 
the rest-frame UV properties of our current sample. In §3, 
we describe the construction of the Lya catalog, discuss the 
method used to measure Lya equivalent widths, and com- 
pute the completeness of Lya detection as a function of lu- 
minosity and redshift. In §4, we discuss the luminosity de- 
pendence of Lya emission in the context of earlier studies 
at z ~ 3. We then turn to the redshift- dependence of Lya 
emission and discuss the important implications of our find- 
ings. Finally, in §5, we examine the extant data on the rate 
of occurrence of Lya emission for candidate sources thought 
to lie beyond z ~ 6. We use this to test the practicality of 
using our test as a valuable probe of cosmic reionisation. We 
summarise the conclusions of our study in §6. 

Throughout the paper, we adopt a A-dominated, 
flat universe with Qa = 0.7, Q m = 0.3 and Ho = 
701i7o kms _1 Mpc _1 . Al l magnitudes in t his paper are 
quoted in the AB system i|Oke fc Gunrj|l983h . 



2 OBSERVATIONS 

We present the results of a new and ongoing Keck spec- 
troscopic survey of photometrically-selected B, V and i! - 
band 'dropouts' in the n orthern and southern GOODS fields 
i|Giavalisco et al.ll2004bl ). The GOODS fields were selected 
for this survey on account of the depth and precision of their 
unique multi-color photometric data useful for selecting tar- 
gets, as well as the availability of associated Spitzer and 



AGN-related properties (Eylcs et al. 


120051; lYan et all 20061; 


Evles et al.l 120071; 


Stark et al. 2007a 


, 2009; Gonzalez et al.l 


20101; Labbe et all 


2009). 



2.1 The Keck/DEIMOS Survey in GOODS-N and 
GOODS-S 

The majority of spectra discussed in this paper were ob- 
tained from an ongoing survey undertaken with the DEep 
Imaging Multi-Object Spectrograph (DE IMOS) at the Nas - 
myth focus of the 10 m Keck II telescope dFaber et al.l l2003). 
DEIMOS is comprised of eight 2k x 4k CCDs spanning 
roughly half of the ACS GOODS field of view (~ Vo'.l x 5'0) 
on the sky. Our first observations have primarily focused on 
targeting the B and V-drop populations. Although future 
DEIMOS observations will expand coverage of higher red- 
shift i' and 2-band dropouts, we include early data taken in 
this territory, some of which is described independently in 
Bunker et al (2010, in preparation). 

The recently-studied target list is primarily selected 
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from the B, V and i' dropout samples discussed in 
IStark et all (|2009h . To this we will add earlier Keck data on 

1 dro ps discussed in §2.2 (jBunker et al 1 l2003l : IStanwav" et al,l 
|2004 Bunker et al. 2010, in preparation) as well as newly- 
discovered 2-drops from the recent WFC3/IR UDF cam- 
paign (see below) . Optical magnitudes range from relatively 
bright systems (zsso — 23.5) to the faintest dropouts ob- 
served in GOODS (zss o ~ 2 7.5). Although the ACS pho- 
tometry in IStark et al.1 (|2009l ) was based upon the GOODS 
version rl.l ACS multi-band source catalogs, for this anal- 
ysis we have updated the photometry to the recently re- 
leased GOODS version 2.0 catalogs. These new catalogs con- 
tain significantly deeper zsso-band imaging. The 5a limiting 
magnitudes vary across the different frames, but are typi- 
cally 28.0 in F435W (B 435 ), 28.2 in F606W (Veoe), 27.9 in 
F775W (1775), and 27.8 in F850LP (z S5 o) when measured in 
0.5 arcsecond diameter apertures. At the median redshifts 
expected for these populations (2 = 3.8 for B-drops, 2 = 5.0 
for V-drops, and z = 5.9 for i'-drops, e.g., Bouwens et al. 
2007), these limits correspond to absolute magnitude limits 
of M uv ~ -18.2, -18.8 , and -19.1 for B, V, and i'-drops. 

A summary of the various observing campaigns is given 
in Table 1. Apart from the i'-drop exploratory survey dis- 
cussed by Bunker et al. 2010, (in preparation), in our 
campaigns during 2004-2007 targets were selected to fill 
empty regions on slitmasks designed for other purposes, 
e.g., studying the kine matics of disk galaxies at 2 ~ 1 
|MacArthur et al]|2008l ). For these observing runs, we used 
the Gold 1200 line mm -1 grating which provided coverage 
between 5570 and 8210 A (allowing Lya to be detected in 
the redshift range 2 = 3.6 to 5.8), with a spectral pixel size 
of 0.3 A pixel -1 . The spectral resolution measured from 
skylines was 1.4 A. 

The main survey began in earnest in 2008, when we 
started a dedicated programme geared at obtaining spec- 
tra of 2 > 3 dropouts. Most spectra were obtained from 
5 slitmasks observed in April 2008 and March 2009 (Ta- 
ble 1). These slitmasks utilised the 600 line mm -1 grat- 
ing on DEIMOS, providing spectroscopic coverage between 
4850 Aand 10150 A (allowing Lya to be detected between 

2 ~ 3.0 and 7.3) with a spectral pixel scale of 0.7 A pixel -1 . 
These spectra provided a resolution (measured from sky- 
lines) of ~ 3.5 A. This setup allowed us to efficently follow- 
up B and V-band dropouts simultaneously. Each individual 
mask was observed for ~ 5 — 7 hours and contained ~ 80- 
110 dropouts. Seeing was typically 0'.'8, ranging between 0'.'5 
and l'.'O. 

We observed a final slitmask toward GOODS-S on 19- 
20 October 2009. For this run, we prioritised i'-drops and 
newly-discovered 2-drops dOesch et al. 2010l; McLure et al.l 
l201oj i iBunker et al.1 [2OI0I ; IWilkins et al.1 l2010al ) and there- 
fore opted for the higher resolution 830 line mm" 1 grating 
blazed at 8640 A with order blocking filter OG550. With 
this set-up, the spectra generally provide coverage between 
6800 and 10100 A (allowing Lya to in principle be detected 
between 2=4.5 and 7.3), with each spectral pixel spanning 
0.46 A. Skylines are measured to have a FWHM of ~ 2.4 A, a 
significant improvement upon the resolution obtained with 
the 600 line grating. However, poor seeing (average of 1.3 
arcseconds) and fog made redshift confirmation particularly 
difficult. 

All data were reduced using the spec2d IDL pipeline de- 



veloped for the DEEP2 Survey. The final reduction provides 
two-dimensional (2D) spectra and variance arrays, along 
with a one-dimensional (ID) extraction at the expected po- 
sition of the dropout. Wavelength calibration was typically 
obtained from Ne+Xe+Cd+HG+Zn reference arc lamps. In 
general, the final wavelength solution is accurate to within 
~ 0.1 A. Examples of reduced ID-spectra are presented in 
Figure [T] 

We flux calibrated the spectra using spectroscopic stan- 
dard stars observed in the 24-26 March 2009 observing run. 
We tested the flux calibration derived from these standards 
using spectra of alignment stars included on the slitmask 
(observed in 2 arcsec x 2 arcsec boxes). We compute opti- 
cal broadband magnitudes for the alignment stars from the 
flux-calibrated spectra using the appropriate filter transmis- 
sion functions. The magnitudes measured from the spectra 
match those from the ACS images to within a factor of ~ 2. 
We bootstrapped a flux calibration on spectra from observ- 
ing runs for which spectroscopic flux standards were not 
taken (03-05 April 2008) using the measured flux in align- 
ment stars that are in common between the 2008 and 2009 
observing runs. 

Using the flux calibrations, we compute our typical la 
flux sensitivity across the DEIMOS spectra for the 600 line 
grating. These measurements predict that we should detect 
continuum at the la level for sources with V — 25.5 (per 
spectral resolution element in ID spectra which have been 
extracted over ~ 1") across much of the wavelength regime 
covered by the 600 line spectra. This prediction is consis- 
tent with expectations based on the optical magnitudes of 
sources that show continuum traces in the spectra. The av- 
erage 10a limiting line flux is 1-I.5xl0 -17 erg cm -2 s -1 be- 
tween 6500A and 9300A, although sky lines become much 
more common toward the red side of the spectra. In a later 
section, we present more detailed simulations that reveal the 
completeness for emission lines of a given flux and redshift. 

For the most recent run (October 2009), the limiting 
sensitivity was significantly worse than in previous runs due 
to the seeing and fog. Based on the signal obtained from 
the bright alignment stars on the mask (each of which has 
known broadband magnitudes) and the signal from our spec- 
troscopic standard stars, we estimate that our continuum la 
sensitivity (computed per resolution element with a spatial 
extraction width of ~ 1.5") was ~ 3.5 — 4. Ox 10 -19 erg cm -2 
s" 1 A -1 in between sky lines (corresponding to a continuum 
magnitude of i' ~ 24.25). Assuming typical line widths, this 
limit translates into a 10cr line flux limit of 2. 5-4. Ox 10 -17 
erg cm -2 s -1 for Lya. As most of the sources targeted in 
this run were faint (2850 > 26 for the i'-drops), this means 
that we are only sensitive to very strong emission lines (the 
lOcr rest-frame equivalent width limit is ~ 75-100 A for 
sources with zgso — 26.5). For fainter sources (e.g. the ma- 
jority of WFC3 2-drops in the UDF), the equivalent width 
limits are too large to enable detection of Lya. 

To summarise the current status of our DEIMOS ob- 
servations, we have obtained 549 spectra of B, V, i! , and 
2-drops. To boost the S/N of our spectra, many dropouts 
were observed on multiple masks; hence, each spectra does 
not correspond to a unique sources. Accounting for this, we 
observed a total of 268 B-drops, 95 V-drops, 19 i'-drops, 
and 17 2-drops. We combine this sample with archival Keck 
and VLT spectra in the following two subsections, and in 
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Number 


Field 


Mask ID 


Date 


t cxp (ksec) 


N B 


N v 


Ni 


N z 


Grating 


I 


ooons-s 


GS031 


08-09 January 2003 


19.8 


o 


o 







1200 


2 


GOODS-N 


GN031 


02-06 April 2003 


37.8 


o 





5 





1200 


3 


GOODS-S 


GS041 


11 December 2004 


21.9 








20 





1200 


i 


GOODS-S 


GS051 


31 October, 01-02 November 2005 8.4 








17 





1200 


5 


GOODS-S 


GS071 


10-11 November 2007 


18.0 


7 


7 








1200 


6 


GOODS-N 


GN081 


03-05 April 2008 


21.6 


85 


9 








600 


7 


GOODS-N 


GN082 


03-05 April 2008 


21.6 


95 


12 








600 


8 


GOODS-N 


GN083 


03-05 April 2008 


20.4 


86 


11 








600 


9 


GOODS-N 


GN094 


24-26 March 2009 


18.0 


45 


63 








600 


10 


GOODS-N 


GN095 


24-26 March 2009 


25.2 


43 


36 








600 


11 


GOODS-S 


GS091 


19 October 2009 


19.2 








19 


17 


830 






Total Uninue Keek/DRIMOS Samnle 




268 


95 


64 


17 






















Total Unique VIT/FORS? Sample 




83 


56 


56 











Comhinerl TTnimip Samnle 




351 


151 


108 


17 

















Table 1. Summary of observations with Keck/DEIMOS and VLT/FORS2. N B , N v , Nj, N z denote the number of B, V, i' , and 2-drops 
observed on each mask. To maximise S/N, many sources were observed on multiple DEIMOS masks; we account for such duplication 
when computing the total Keck sample size. The i'-drops in the first two rows were originally published in Bunker et al. (2003) and 
Stanway et al. (2004), and the i'-drops presented in rows 3 and 4 are from Bunker et al. (2010, in preparation). The VLT FORS2 sample 
is take from observations presented in Vanzella et al. (2009). The combined total sample in the last row includes all unique dropouts in 
the Keck and VLT surveys. 



§2.4, we present the absolute magnitude distribution of the 
combined VLT and Keck spectroscopic samples. 



2.2 Archival Keck Spectroscopy in GOODS-N/S 

Additional i'-drops were observed with Keck/DEIMOS be- 
tween 2003 and 2005 (Table 1). Early observations were pre- 
sented by Bunker et al. (2003) and Stanway et al (2004), and 
an updated discussion, including additional DEIMOS obser- 
vations from 2003-2005, is given by Bunker et al. (2010). 
These sources were observed with the 1200 line mm - (de- 
scribed above) with seeing of 0'.'7 - l'.'O. 

In total, 45 i'-drops were observed in the two GOODS 
fields on 4 separate slitmask s. Of the 45 sour ces observed, 
12 were included in the IVanzella et all (|2009h VLT/FORS 
observations discussed below. Integration times ranged be- 
tween 2.3 and 10.5 hrs. The Z85o-band magnitudes of the 
sources ranged between 24.7 and 28.3. Given the consider- 
able range of exposure times, we take care to estimate the 
equivalent width completeness for sources of different mag- 
nitudes on the various masks that were observed. 



2.3 Archival VLT/FORS Spectroscopy in 
GOODS-S 

Two programmes aimed at following up high-redshift 
dropouts in GOODS-S have been conducted with the VLT. 
The first of these used the FORS2 multi-object spectrograph 
(Vanzella et al. 2005, 2006, 2008, 2009), a nd the second 
used t he VIMOS multi-object spectrograph iBalestra et al.l 
(2010). Both teams have released their datasets to the pub- 
lic. In this paper, we focus on the FORS2 survey, as its sur- 
vey characteristics (resolution, spectral coverage) are closest 
to the DEIMOS survey. We discuss the basic FORS2 survey 
details below. 

Between September 2002 and October 2006, the VLT 
FORS2 multi-object spectrograph was used to observe 



sources identified in the GOODS imaging of CDF-S. In to- 
tal, 38 FORS2 masks were obtained using the 3001 grism 
without an order-separating filter. Each mask was observed 
for roughly 4-6 hours. In general, the spectra provide cover- 
age between 6000 A and 10000 A, with a spectral resolving 
power of R=660 which provides resolution of ~ 13 A at 
8600 A. 

The FORS2 database provides redshift classifications 
and quality grades (ranging from A to C) for the en- 
tire spectroscopic sample. The colour criteria used to se- 
lect the dropou t s in th e FORS2 sample are discussed in 
iGiavalisco et alJ (|2004al ) : however, for the B-drops, a slight 
variation in the colo ur-selection was adopted (for details see 
IVanzella et"alll2009h . 

In general, these criteria are very similar to those we 
have adopted for our DEIMOS survey. Using the coordi- 
nates provided in the public FORS2 database, we query the 
version 2.0 ACS catalogs for GOODS-S and measure opti- 
cal magnitudes in an identical fashion as described in §2.1. 
Adopting the selection criteria used for the DEIMOS sam- 
ple (see Stark et al. 2009), we find 83 B-drops, 56 V-drops, 
and 56 i'-drops that were observed with FORS2. Redshifts 
were obtained for 48 B-drops (46 with z > 3), 37 V-drops 
(32 with z > 4), and 26 i'-drops (21 with z > 5). 

The magnitude distribution of the FORS2 sample is 
generally weighted toward brighter sources, with few B and 
V-drops with magnitudes fainter than 2850 — 26. Given the 
inherent faintness of i'-drop samples, the magnitude distri- 
bution is significantly fainter than for the lower redshift 
dropout samples, with zsso — 25 — 27. In Figure we 
plot a comparison of the absolute magnitude distribution 
of the FORS2 and DEIMOS dropout samples. The current 
DEIMOS sample contributes 76% of the B-drops and 63% 
of the V-drops and crucially extends the spectroscopic cov- 
erage to lower luminosities. 
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Figure 1. Montage of DEIMOS ID spectra from our survey arranged according to redshift and luminosity. B-drops are shown in the left 
column, and V-drops appear in the right column. Each panel contains a label indicating the apparent magnitude of the continuum and 
redshift. Spectra are smoothed to a spectral pixel size of ~ 3 A. Where present, dashed vertical lines denote Lya, Si II 1260, OI+SiII 
1303, and CII 1334 lines. 



2.4 Final Spectroscopic Sample 

When combined, the FORS2 and DEIMOS surveys contain 
a total of 627 spectra of unique high-redshift dropouts (351 
B-drops, 151 V-drops, 108 i'-drops, and 17 z-drops). In Fig- 
ure O we present colour-colour diagrams of the entire B 



and V-drop spectroscopic sample, indicating the selection 
criteria. The i'-drops were selected using the standard sin- 
gle i' ~ z > 1.3 colour criterion (see Figure [3). Comparing 
the distribution of selected ta rgets on these figur es with the 
larger sample of sources in the lStark et alj ([2009) photomet- 
ric catalog, we demonstrate that the colours of our targets 
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Figure 3. Colour-colour diagrams for the B-drops (left) and V-drops in combined VLT and Keck spectroscopic surveys. Objects with 
Lya in emission are denoted with a second red circle. The dotted lines show the Lyman break selection cr i teria adopted in this paper. 
The small grey circles show the distribution of colours for the parent population of LBGs from St ark et al 
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Figure 2. Absolute magnitude distribution of FORS and 
DEIMOS spectroscopic samples. The new DEIMOS B, V, and 
-drops are shown in shaded blue, and below the FORS2 observa- 



tions jVanzella et al]|2009t) are shown in light red. The DEIMOS 
datasets comprise the majority of the spectra considered in our 
analysis. 




Figure 4. i' — z colour versus magnitude for i'-drops with spec- 
troscopic observations with the dashed line showing the selec- 
tion cut. Points surrounded by a red open circle indicates sources 
with Lya seen in emission, and the small grey circles indicate the 
colours and magni tudes of i'-drops from the parent LBG samples 
l|Stark et alJl200Sh . 



with spectroscopic coverage span nearly the full range of the 
parent LBG population. 

Figure [5] also shows the apparent magnitude distribu- 
tions of the various dropout samples. While previous surveys 
have focused on sources with apparent magnitudes brighter 
than m ~ 26, a major achievement of our survey is that we 
have been able to push significantly below this limit with 
DEIMOS. Our primary motivation for doing so is that the 
characteristic UV luminosity shifts to lower values at higher 
redshifts (e.g., Bouwens et al 2007), so to compare source 
properties over 3 < z < 7 in a meaningful manner, probing 
deep becomes a necessity. To demonstrate our survey has 
achieved this over 3 < z < 6, in Figure [6] we compare the 
absolute magnitude distribution of our various dropout sam- 
ples. Current estimates of the chara cteristic absolute mag- 
nitude at z ~ 7 are M uv ~ -19-9 (|Bouwens et al.l l2009al : 
lOuchi et ai1l2009l ; [Oe"sch et alj|2010h . As objects more lumi- 



nous than this are exceedingly rare at z > 6 — 7, if we are to 
make consistent Lya fraction comparisons, we must probe 
to at least this depth over 4 < z < 6. 

Redshifts were determined via visual examination of the 
spectra (see Figure [7] for a redshift distribution). For UV 
faint sources, the continuum is too faint for identification of 
the Lyman break or absorption lines, so we can only measure 
spectroscopic redshifts for those sources with Lya in emis- 
sion. For UV bright sources, redshift identification is per- 
formed via a combination of the Lyman break, interstellar 
absorption lines, and/or Lya emission. We classify all red- 
shifts according to their quality, ranging from A (definite), 
B (secure), C (possible/likely). Nearly all Lya emitters fall 
into the first two categories owing to the combination of line 
profile, lack of other emission lines, and strong continuum 
break (in spectra and/or imaging), but some of the absorp- 
tion line detections are much more tenuous owing to sky 
line residuals and noise. As the results of this paper are in- 
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Figure 6. Distribution of absolute magnitudes of sources in the 
combined DEIMOS and FORS2 spectroscopic survey (shaded 
blue). The shaded light red, green, and grey denote the mag- 
nitude distribution of the B, V, and i'-drop samples. 




Figure 7. Rcdshift distribution of B, V, and i'-drops with spec- 
troscopic confirmation in GOODS-N and GOODS-S (dark blue). 
This sample includes absorption line systems in addition to the 
Lya emitters discussed in this paper. The total sample contains 
179 galaxies with rcdshifts in the range 3 < z < 4.5, 87 galaxies 
with redshifts over 4.5 < z < 5.5, and 44 galaxies with rcdshifts 
over 5.5 < z < 6.5. For reference, we display the redshi ft distri- 
bution of the VLT/FORS2 survey l lVanzella et alj|2009h in light 
red. 



dependent of the absorption line sample, we delay further 
discussion of the absorption line catalog until a subsequent 
paper (Stark et al. 2010, in preparation), while details of the 
Lya selection are described in the following section. How- 
ever, we note that, in total, the FORS2+DEIMOS sample 
contains 179 redshifts over 3 < z < 4.5, 87 redshifts over 
4.5 < z < 5.5, and 44 redshifts over 5.5 < z < 6.5. 



3 EMISSION LINE MEASUREMENTS 

3.1 Constructing the Lya catalog 

Emission lines were initially identified visually in the 
two-dimensional DEIMOS spectra and later in the one- 
dimension extractions at the position of the dropouts. We 
took care to distinguish Lya from other emission features 
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Figure 8. Spectra of a J125 — 27.3 z-drop showing tentative Lya 
emission at 9044 A. Top panel: The 2D spectrum shows a sig- 
nificant emission feature detected at the la level (black corre- 
sponds to positive flux in the image) in between two sky lines 
and cen tered exactly at the position of z-drop '3' in lWilkins et al] 
l|2010ah . Bottom panel: The Id spectrum of the galaxy and asso- 
ciated noise. The spectrum and noise are not smoothed to avoid 
blending with skylines at 9038 A and 9048 A. The emission 
line between these features lies in a region of low noise, spanning 
~ 10 A in width. 



(e.g. [Oil] which is resolved by DEIMOS) which correspond 
to lower-redshift galaxies. Lya emission is detected in 152 
of the DEIMOS dropouts. 

For the VLT/FORS data set, the redshifts were deter- 
mined in IVanzella et all lj2009h . As discussed in §2.3, to en- 
sure a uniform selection across the Keck and VLT samples, 
we performed our own photometric selection. With the re- 
sulting subset of 193 galaxies, we identified those objects 
showing Lya in emission in their one-dimensional spectra. 

We defined Lya redshifts (zL ya ) for each object as the 
wavelength at which the Lya line is at its peak flux value. 
As Lya is generally redshifted by at least ~ 300 km s -1 
from the frame of rest of the galaxy (e.g. IShaplev et al.l 
120031 ; IVanzella et ail 120091 ). we note that this redshift is 
not necessarily equivalent to the systemic redshift of the 
galaxy. Across both surveys, we identified Lya emission 
in 108 B-drops, 63 V-drops, and 28 i'-drops. The relative 
colour distribution of these sources in relation to the larger 
spectroscopically-targetted sample is shown in Figures[3]and 

m 

In addition to the B, V, and i'-drops, we present tenta- 
tive spectroscopic confirmation for one of the z-drops that 
was studied in poor conditions in the October 2009 DEIMOS 
run (Table 1). The only z-drop to show a clear emission 
feature at the expec ted object location was source 3 in 
IWilkins et all l|2010al ). This object, identified in the GOODS 
ERS WFC3 imaging of CDF-S (PI: O'Connell), is reason- 
ably faint (J125 — 27.3) but is brighter than most of the 
z >, 6 sources detected in the UDF (e .g. lOesch et al.ll20ldl : 
iBunker et al.ll20ld ; lMcLure et al.ll2010T ). The DEIMOS spec- 
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Figur e 9. Spectral energy distribution of z-drop '3' from 
IWilkins et al.l l|201 0a*l. The datapoints show detections and up- 
per limits from ACS and WFC3 imaging of CDF-S. Overplotted 
is a Bruzual-Charlot population synthesis model with a redshift 
fixed at z = 6.44, the spectroscopic redshift inferred from the 
emission line detected in Figure [8] The broadband imaging data 
show a prominent break where the Lyman break is predicted and 
a significant excess of flux in the medium-band Y098 filter (rela- 
tive to the Ji25-band), as expected in the presence of strong Lya 
emission. 



trum illustrates a significant 7a emission line in between 
sky lines and centered at 9044 A (Figure [8]). The S/N and 
nearby sky line do not enable a robust measure of the asym- 
metry of the line. The broadband SED is well fit by a source 
z ~ 6.44 with strong Lya contaminating the medium-band 
y"o98 filter (Figure [9]). If this emission feature is Lya, it would 
lie at z — 6.44, consistent with the photometric redshift pre- 
dictions. Our flux calibration suggests a flux of 1.5 x 10~ 17 
erg cm" 2 s" 1 and hence an equivalent width (rest-frame) of 
83 A. We will confirm the redshift of this source in future 
campaigns of CDF-S. 



3.2 Computation of Lya properties 

Next we compute the flux (FL ya ) and equivalent width 
(Wljo) for each Lya emitter in our spectroscopic sample. 
Previous attempts to constrain evolution in the prevalence 
of Lya emitters have focused on measuring evolution in the 
luminosity function of Lya emitters (selected via narrow- 
band filters). For our Lya fraction test described in §1, we 
are interested in determining the percentage of LBGs of a 
given luminosity with Lya emission much stronger than the 
continuum flux. This test thus relies on accurate measure- 
ments of the Lya equivalent widths (W Lya)- 

In principle, Wlj/q provides a more robust constraint 
than FLya as it does not rely on an absolute flux calibra- 
tion. However since the majority of spectra in our sample 
have broadband continuum magnitudes fainter than the la- 
continuum flux limit of our survey (m ~ 25 — 26, §2.1), we 
can only place an upper limit on the continuum flux mea- 
sured in this subset of our spectra, thus producing a lower 
limit to the observed WLya- We can obtain a better es- 
timate of the equivalent width by adopting the continuum 
flux measured from broadband imaging, taking care to avoid 
filters that are contaminated by Lya. We detail the specific 
procedure used to measure Wi ja below. 
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For each spectrum, we compute the line flux, Ftya., by 
summing emission in excess of the continuum between 1213 
and 1221 A, taking care to avoid contribution from nearby 
sky lines or artifacts and ensuring that the spatial extraction 
box covers the entire spatial width of the line emission. 

The measurement of the continuum flux depends on 
whether or not continuum is detected. For those objects 
with continuum detections, we compute the average contin- 
uum level redward of Lya, c re[ j, by averaging the flux in re- 
gions between OH sky lines between 1225 A a nd 1255 A in 
the re st-frame of the galaxy. As discussed in iKornei et al.l 
(2009), this should minimize contribution from nearby ab- 
sorption features. The observed-frame Wi S o is then com- 
puted by taking the ratio of F^ya and c re d- 

For those objects without continuum detections, we 
compute equivalent widths using the continuum level just 
redward of Lya with the broadband photometry discussed 
in the previous section and taking care not to include the 
contribution of Lya to the continuum. For the B-drops and 
V-drops, we use the 1775 and zsso-band fluxes, respectively. 
These measurements provide the luminosity at rest-frame 
~ 1500 A, rather than just redward of Lya. The vast ma- 
jority of sources are very blue (and hence nearly flat in f„), 
requiring no correction between the flux at 1500 A and 
that at 1240 A (the central value in the wavelength-region 
summed for continuum measurements above). For those few 
sources with very red UV slopes (measured from the broad- 
band SEDs presented in Stark et al. 2009), we apply a 
small adjustment to correct for the change in flux between 
1500 A and 1240 k. 

For galaxies with redshifts above z ~ 5.7, the zsso-band 
filter is contaminated by both Lya emission, if emission is 
present, as well as Lya forest absorption shortward of rest 
1216 A. Accurate continuum determinations must account 
for this. If sources contain Y or J-band detections via WFC3 
or ISAAC, then we use these measures for the UV continuum 
However, at the moment, the majority of UV-faint i'-drops 
are not detected with significance in the near-IR, so we must 
use zg5o detections for these sources. For galaxies with no de- 
tected Lya emission, we must only correct for the Lya forest 
absorption, which we do assuming the redshift is identical 
to its value inferred from BPZ and that its SED is flat in f„. 
Corrections range from ~ 0.2 mag at z = 6 to 0.7 mag at 
z ~ 6.4. But for z > 5.7 objects with strong Lya emission, 
we also subtract off the line contribution when computing 
the continuum flux. Typical line contributions are ~ 0.1-0.3 
mag. 

To test the reliability of the above method, we examine 
the subset of galaxies with bright UV-continua, comparing 
the continuum flux determined from broadband-imaging to 
that extracted from the spectra. We find that the median 
fractional error using the estimate from broadband imaging 
is ~ 20% (Figurc fTIJ)) . We add this error in quadrature to the 
photometric error on each equivalent width measurement. 

3.3 Completeness of Lya detection 

In order to properly assess the fraction of Lya emitters in our 
sample as a function of UV luminosity and redshift, we must 
consider how completeness varies with apparent magnitude 
and wavelength of the Lya emission line. We estimate the 
typical completeness by adding fake Lya emission to random 
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Figure 10. Reliability of WL ya measurements. The x-axis shows 
equivalent widths (observed-frame) determined directly from 
spectra for sources bright enough for a continuum measurement. 
The y-axis shows the equivalent widths measured for the same 
sources using the continuum flux from the broadband photome- 
try. The dashed lines demarcate values 0.5 X and 2x the WLya 
measured directly from the spectrum. The median fractional er- 
ror from using equivalent widths derived using the continuum flux 
from broadband imaging is ~ 20%. 

positions across the DEIMOS and FORS2 spectra. In each 
case, we measure the line properties (WLya, ^L ya , S/N) of 
the fake emission feature. We ran enough trials (using a 
large number of spectra from each mask) to obtain reliable 
estimates of the completeness of the Lya recovery. 

Our goal is to identify a lower flux threshold above 
which Lya is highly complete for most objects in our sample 
in order to minimize the necessary completeness corrections. 
We illustrate the results of the completeness simulations for 
the 600 line mm" 1 grating on DEIMOS in Figure [TT] The 
masks on this grating included B and V-drops which span 
the redshift range 3.5 < z < 5.5. The simulations reveal that 
sources with hiab — 27.0 and with rest-frame Lya equiva- 
lent widths (Whyafi) in excess of 50 A are generally re- 
covered over nearly the entire redshift range probed by our 
dropouts (see Figure Hip . The high completeness for very 
strong Lya emitters amongst bright sources (m < 26) arises 
as these bright emission features correspond to fluxes that 
are in excess of 10-20<r and are recovered even when they lie 
on top of sky lines. The Whyafi > 50 A emission line com- 
pleteness begins to decline for very faint sources (i' ~ 27.5) 
at the high-z tail of the V-drop redshift distribution. For the 
FORS2 data, the spectra show similarly high completeness 
for B and V-drops with WLya,o > 50 A. These simulations 
suggest that the measured Lya fraction of sources in the 
redshift range 3.5 < z < 5.5 should not suffer from signifi- 
cant incompleteness for B and V-drops more luminous than 
Muv ^ -19. 

For the i'-drops, the completeness in the deep DEIMOS 
and FORS2 spectra remains high for strong (e.g. WL ya ,o > 
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Figure 11. Emission line recovery completeness. The complete- 
ness is determined via simulations placing fake lines in spectra 
and testing the rate at which sources about the equivalent width 
threshold are recovered. From top to bottom, the figure shows 
completeness versus redshift for WL yQ ,o sources with continuum 
magnitudes (longward of the Lyman break) of 26.0, 26.5, 27.0, 
and 27.5. 

50 A ) Lya lines. However, it is clear from Figure [TT] that in- 
completeness is not negligible, particularly for faint (zsso > 
27.0) sources at z > 5.8 which are less than 80% complete. 
We thus adopt this as our z'-drop magnitude threshold for 
inclusion into the Lya fraction test, limiting us to sources 
more luminous than Mtjv — —19.7. We also do not include 
the 17 i'-drops on the mask GS051 (see Table 1) from Bunker 
et al. 2010 (in preparation) owing to its reduced integration 
time of 2.3 hours. 



3.4 Contamination in spectroscopic samples 

Although foreground emission line sources (e.g. [O II] emit- 
ters can be readily distinguished from those revealing Lya 
(§3.1), our Lya fraction test requires that we have a reliable 
sample of LBGs for which no emission line is seen. 

While the Lyman break selection criteria are chosen 
to minimize the inclusion of low-redshift and stellar con- 
taminants, it is clear that interlopers still populate dropout 
samples. For our purposes, it is important that we consider 
the luminosity and redshift dependence of contaminants, as 
these could create artificial trends in our derived Lycv frac- 
tions. We expect \aw-z contamination to increase toward the 
fainter end of our sample. This is primarily due to photo- 
metric scatter; since faint sources are detected with lower 
S/N and have less dynamic range available to constrain the 
break, it is more likely that a faint low-z source will be scat- 
tered into the LBG selection window. However, if the typ- 
ical contamination is fairly low in all luminosity bins (e.g. 
< 10%), then this would only require a minor correction to 
the luminosity dependence of the Lycv fraction. 

To investigate this in more detail, we compute pho- 
tometric redshift probability distributions for the spectro- 
scopic sampl e using the o bserv ed photometric catalogues 
compiled in IStark et al.l (|2009l ). updated to include v2 
GOODS ACS photometry, and the BPZ software (Benitez 
et al. 2000). Deta ils of the photome tric redshift methodol- 
ogy is discussed in lStark et al.l |2009l 'l. Using the probability 
distributions derived from BPZ, we compute the probability 



that each object lies outside the redshift range constrained 
by our spectroscopic observations (typically z < 3.4). We 
then place those galaxies without spectroscopic redshifts in 
bins of UV luminosity and compute the total contamina- 
tion fraction as a function of UV luminosity. As expected, 
the results show that the contamination fraction increases 
toward lower luminosities. We find negligible contamina- 
tion for bright sources (—22 < Mtjv < —20) with SEDs con- 
strained by high S/N photometry. In the two faintest abso- 
lute magnitude bins considered (Mtjv = — 19 and -18), the 
predicted contamination increases to ~ 10%. This suggests 
that, owing to low- 2: contamination, the true Lya fractions 
in the two faintest bins should be 1.11 x greater than de- 
rived. We consider the effects of this in our discussion of the 
Lya fraction in §4.1. 

We now examine whether the SEDs predict contam- 
ination should vary strongly with redshift. As above, we 
measure the contamination fraction (of objects that are 
not Lya emitters) implied for B, V, and z'-drop samples 
from the photometric probability distributions. As above, 
we find negligible contamination at the brighter magni- 
tudes. Combining sources with UV luminosities spanning 
—20 < Mtjv < —18, we find that the contamination fraction 
increases from 2% for V and i'-drops to 5% for B-drops. This 
suggests that in this luminosity regime, the B-drop Lya frac- 
tion will be underestimated by a factor of 1.05, while the V 
and i'-drop Lya fraction will be underestimated by a factor 
of 1.02. Hence low-z contamination will cause the positive 
evolution in the Lya fraction with redshift to be overesti- 
mated by a factor of 1.03 in this luminosity range. As we 
will show in §4.2, while the redshift-dependent trends in the 
Lya fraction are small in their amplitude, this contamina- 
tion effect contributes little to the observed variation with 
redshift. 



4 ANALYSIS 

We now turn to the key questions posed in §1. Armed with 
a large sample of LBGs for which some fraction, XL ya , show 
Lya emission, we discuss what can be learned about the de- 
mographics of line emission in the star-forming population 
and how such trends might affect our ability to use Lya 
as a tracer of reionization. In this analysis section, first we 
present the results, both in terms of the fraction of LBGs 
showing line emission as a function of luminosity and red- 
shift. We then discuss the physical factors that might explain 
these trends prior to the use of Lya emission statistics as a 
possible probe of reionisation. 



4.1 The Luminosity-Dependence of Lya emission 
at high-redshift 

First, we discuss the relationship between Lya emission 
and luminosity in our spectroscopic sample of Lya emit- 
ters between 3 < z < 6. If the Lya fraction varies strongly 
with luminosity, as may be expected given recent claims 
of luminosity-dependen t dust obscuration at h igh-redshift 
l|Reddv fc Steidell 120091 ; iBouwens et al. 2009b). then care 
must be taken to compare only galaxies of similar luminos- 
ity when searching for evolution in the Lya fraction with 
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redshift. As our survey probes to considerably lower lu- 
minosities than past spectroscopic LBG samples at z > 3 
(Muv ~ —18, Figure [6}, our sample is well-suited to inves- 
tigating such a relationship. 

There have been a number of previous studies which ex- 
amine how Lya line strength varies with luminosity. Many 
of these studies have reported a cor relation between Lya 
equivalent width and UV luminosity. IShaplev et al.l (|2003l ) 
binned their sample of z ~ 3 galaxies with Whyafi > 20 A in 
three groups of apparent UV luminosity and found that the 
mean of the Whyafi distribution increased toward fainter UV 
luminosity. Others have examined Lya equivalent widths as 
a function of UV continuum luminosity, revealing a deficit 
in large equivalent width Lya lines in the most luminous 
continuum sources ( Ando et al.l l200d; lOuchi et all 120081; 
Pentericci et all 12003 ; IVanzella et all 120091 ; iBalestra et all 



20101 ). Additionally, by comparing the UV luminosity func- 
tion of LBGs w ith that of narrowband-selected LAEs, 
lOuchi et alj (|2008l ) has shown that Lya emitters are likely to 
be more prevalent at the faint-end of the luminosity func- 
tion. As we discussed earlier, these results are consistent 
with simple theoretical expectations in which low luminos- 
ity galaxies are less obscured by dust (due to lower metal- 
licities) and perhaps have lower column densities (or cover- 
ing fraction) of HI surrounding them (e. g. IVerhamme et all 
l200d . l200Sl : ISchaerer fc Verhammell200Sl ). Perhaps of equal 
or greater importance is the bulk veloc ity field of the HI 
jShaplev et al J 120031 ; Isteidel et al.ll2010h . which we discuss 
further in §4.3. 

However, others have found the evidence for a corre- 
lation between lumi nosity and Lya equivalent width to b e 
less convincin g Ce.g.lSteidel et alJhoOOl : iNilsson et al.ll2009l ). 
In particular. INilsson et al. l (|2009l ) "demonstrated that the 
dearth of luminous galaxies with extreme Lya emission 
in magnitude or flux-limited surveys does not require the 
equivalent width distribution to be luminosity-dependent. 
The lack of such extreme emitters is actually expected in 
a magnitude or flux limited survey, as the most luminous 
and the most extreme emitters are both rare, causing this 
portion of parameter space to poorly represented unless 
very large volumes are covered. This is an important re- 
alisation but only considers LBGs with continuum lumi- 
nosities brighter than Muv — —20, significantly more lu- 
mino us than the feeble s ources probed in our survey. Like- 
wise, iKornei et al.l (2009) found only marginal evidence for 
a correlation between Whya and UV continuum luminosity 
in a large sample of z ~ 3 LBGs, but concluded that this 
may result from the limited dynamic range in UV luminosity 
probed by their sample. 

Figure [12] shows the distribution of rest-frame equiva- 
lent width, Wl hq ,o for Lya emitters in our survey as a func- 
tion of rest-frame UV luminosity, Muv- In considering this 
figure we must choose a sufficiently bright Wlj 0| o threshold 
to avoid incompleteness in the faintest sources. As demon- 
strated in §3.3 and Figure [TT1 sources with Wljo.o in excess 
of 50 A are detected with high completeness (> 90%) across 
our spectra, so we adopt this value as our equivalent width 
threshold. The data show an apparent lack of strong line 
emission among t he most luminous dropouts, as has been 
found elsewhere (lAndo et al ] 120061 : IPentericci et al.1 120091: 
IVanzella et"ai] f2009). yet as demonstrated in INilsson et al.l 
( 20091 ). this does not necessarily imply that Lya emission is 
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Figure 12. Rest-frame Lyo equivalent width (WL ya ,o) as a func- 
tion of rest-frame absolute magnitude, Muv f° r Lya emitters in 
our spectroscopic sample. The sample is nearly complete at all 
Muv f°r sources with WL ya ,o > 50 A (dotted line). 



less common in luminous galaxies. However, when Figure IT2"1 
is viewed in concert with the absolute magnitude histogram 
of our spectroscopic sample (Figure [6}, it becomes clear that 
strong line emission must be more common in low lumi- 
nosity dropouts, for the number of low luminosity galaxies 
(Muv > —19.5) targeted in our campaign is as low as the 
most luminous sources (Muv < —21.5), but the number of 
strong line emitters is far greater among the feeble sources. 
We emphasize that owing to our equivalent width threshold 
of Wi Ba ,[) ~ 50 A, the Lya emitters amongst low luminosity 
sources are secure and robust (S/N > 10) detections (simi- 
lar to those at the top of Figure 1); hence this result is not 
driven by noisy features. 

In order to most clearly quantify the luminosity- 
dependence of Lya emission, we must compute the frac- 
tion of LBGs in our spectroscopic sample that show strong 
Lya emission, Xhya, as a function of emerging UV luminos- 
ity. Given that the majority of our spectroscopic sample is 
fainter than our continuum flux sensitivity (Figure [SJ, it is 
not possible to measure redshifts for faint sources without 
emission. Hence if we were to compute XL yQ based soley on 
sources with confirmed redshifts (as has been done previ- 
ously for brighter samples), we would artificially increase 
%Lya toward lower luminosities. To avoid this bias, we de- 
fine XL ya as the number of Lya emitters above some Wi Ba ,o 
threshold divided by the total number of dropouts placed on 
our slitmasks. 

The error on iLya is computed as follows. We first de- 
rive the Poisson error from the number of sources considered 
in each luminosity or redshift bin. In addition to the Pois- 
son error, each measurement is subject to additional error 
owing to uncertainty in the equivalent width measurements. 
We thus conduct Monte Carlo simulations, randomly vary- 
ing the equivalent width of each galaxy assuming a normal 
distribution with mean and standard deviation correspond- 
ing to the measured WL ya ,o and aw- For each realization, 
we compute the luminosity-dependent Lya fraction. Con- 
sidering all the trials, we compute the standard deviation 
in the XL ya distribution. We then combine this error term 
in quadrature with the random error derived above. In do- 
ing these simulations, we also consider whether equivalent 
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width error may artifically scatter a net flux of feeble galax- 
ies above our equivalent width threshold, introducing the 
observed trend. However the simulations demonstrate that 
this is not the the trend is readily apparent in each 

realisation. 

In constructing the Lya fraction, we must also recog- 
nize that the wavelength coverage of some spectra is such 
that Lya would not be recovered across the full redshift 
range over which those dropouts might be selected. This 
applies mostly to FORS2 spectra of B-drops which can only 
detect Lya for sources with z > 3.9. Given that B-drops 
are expected to span the redshift range 3.5 < z < 4.5, this 
suggests that a large number of B-drops would not be re- 
covered in the FORS2 spectra even if they showed Lya in 
emission. We therefore limit our study of B-drops to those 
with spectra that allow Lya to be detected over the entire 
redshift range expected. This is a particular advantage of 
the DEIMOS component of our survey where the numer- 
ous B-drops taken with the 600 line mm -1 grating are fully 
sampled. 

Following this procedure, we compute the luminosity- 
dependence of our entire spectroscopic sample. We find that 
iLya is considerably larger in low luminosity galaxies (Figure 
I13p . increasing from 10-20% for luminous sources (Muv — 
-21) to 60-70% for feeble galaxies (M uv ~ -18-19). A sim- 
ilar trend is seen if we adopt larger equivalent width thresh- 
olds. The Lya fraction rises slowly with decreasing luminos- 
ity over —22.0 < Muv < —20.5 but then begins to increase 
more rapidly at lower luminosities( —20.5 < Muv < —18.5). 
We thus fit first order polynomials over each of these lumi- 
nosity ranges, finding iLyo = 1-09 + 0.047 Muv at the lu- 
minous end and iLya = 5.46 + 0.26 Muv for the lower lumi- 
nosity sources. Given that the Lya fraction increases slowly 
for luminous sources, it is perhaps no surprise that previous 
studies which were limited to this luminosity regime found 
only moderate evidence for luminost y dependent trends 
l|Nilsson et al.ll2009l ; iKornei et all|2009j ); it is only by prob- 
ing lower luminosity galaxies that we start to see a clear 
trend in the prevalence of Lya with luminosity. 

Intriguingly, the "break" in the luminosity-dependence 
of the Lya fraction in Figure [13] occurs at very similar lumi- 
nosities to the characteristic UV luminosity o f the rest-frame 
UV 4 < z < 6 luminosity function (e.g., iBouwens et al.l 
120071 ). At luminosities greater than LJ} V , the Lya fraction 
is low and increases slowly, but below this luminosity, the 
Lya fraction increases to much larger values. While perhaps 
coincidence, this may suggest that whatever process modu- 
lates the knee of the luminosity function at z > 4 may play 
a role in the escape of Lya photons. 

In order to put these results in context, it is interest- 
ing to estimate the escape fraction of Lya photons that is 
implied by these large equivalent widths. In particular, we 
are interested in the Lya escape fraction that is implied by 
the WLya,o — 50 A threshold we have adopted. Assuming 
a Salpeter IMF with M uppGr =120 Mq, case B recombina- 
tion, constant star formation, and metallicity ranging be- 
tween Z=l/20 Zq and Zq, Lya equivalent widths can be 
as high as ~ 200 — 300 A in the first few Myr of the star 
formation episode, asymptoting to ~ 10 A after 10 Myr 
l|Malhotra fc Rhoadsll2002i ; ISchaererll2003h . Hence our equiv- 
alent width threshold corresponds to a Lya escape fraction 
of 15-50%. Of course, this assumes that all ionising photons 



are absorbed; if there is significant Lyman continuum leak- 
age, then the predicted equivalent widths would decrease, 
increasing the inferred Lya escape fraction. Alternatively, if 
dust is confined to cold neutral clouds, then the maximum 
Lya equivalent widths may be larger than quoted abov e 
(|Neufeldlll99ll : [h anscn k Oh 2006; iFinkelstein et al.ll2008r i. 
decreasing the implied Lya escape fraction. 

These results provide clear evidence that Lya emission 
becomes continuously more prevalent among lower luminos- 
ity star forming galaxies. Indeed, it appears that the major- 
ity of feeble (Muv < —19) sources are strong Lya emitters. 
These results appear to indicate that the escape fraction of 
Lya photons (relative to that of far-UV continuum photons) 
is strongly luminosity-dependent and at very low continuum 
luminosities may commonly exceed the 5% robustly derive d 
(at z ~ 2) via Lya and Ha surveys (|Haves et al.ll2010h . 
In the following section, we examine whether similar trends 
are seen with redshift, while in §4.3 we examine the physi- 
cal mechanism is that are likely to be governing the escape 
fraction of Lya photons at z > 3. 

4.2 Variation in the Lya fraction in the redshift 
range 3 < z < 6 

We now examine the redshift evolution of the prevalence of 
Lya emitters in the Lyman break galaxy population over 
3 < z < 6. Since the IGM appears to be highly ionised 
over this redshift range, this measurement provides the op- 
portunity to understand the extent to which factors other 
than IGM attenuation affect the Lya fraction. By calibrat- 
ing these effects over this redshift interval, we can more ac- 
curately detect the signal of reionisation on the Lya fraction. 

We consider luminosity-dependent samples in two sep- 
arate redshift bins, 3.0 < z < 4.5 and 4.5 < z < 6.0 (Fig- 
ure [13}. Sources without spectroscopic redshifts are placed 
into one bin or the other based on their photometric red- 
shift (see lStark et al.|[2009l for discussion of photometric red- 
shifts). In each UV luminosity bin, the fraction XLya in- 
creases with redshift. The two bins with the lowest error 
and incompleteness (Muv — —21 and -20) show increases 
of 60% and 40%, respectively. Our equivalent width thresh- 
old and redshift binning ensures that this result is not biased 
by redshift-dependent incompleteness. To determine the dif- 
ferential growth, we compute the average change in Lya 
fraction, AxL ya , across all luminosity bins (weighted by the 
inverse variance of each bin) and compute Az using the me- 
dian redshift in each of the two bins. With this approach, we 
find that the Lya fraction increases with redshift following 
dxLy a /dz ~ 0.05. As we mentioned in §3.4, low-z contami- 
nation appears to decrease very slightly with redshift. How- 
ever, we find that this effect introduces such small changes in 
the Lya fractions (< 1%) in each redshift bin, such that the 
weighted Lya fraction redshift evolution remains as quoted 
above. Hence Lya fractions grow by nearly Ax-L ya ~ 0.1 at 
fixed Muv between z ~ 4 and z ~ 6. 

We emphasize that this trend is not driven by biases 
associated with LBG selection. It is well-established that 
the presence of strong Lya can affect the broadband c olours 
used for dropout selection |Stanwav et al.ll2007l . [2008T) . Line 
emission can either boost the dropout color (if the redshift 
is at the high end of the distribution with Lya in the red- 
der filter) or it can dilute the color (if the redshift is at 
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Figure 13. (Left:) Fraction of spectroscopic dropout sample showing strong Lya emission (WL ya ,o > 50 A) as a function of UV 
luminosity. The dashed lines corresponds to first order polynomial fits to the Lya fractions in the range —22.0 < A'/tjv < —20.5 and 
—20.0 < Mtjv < —18.5. (Right:) Fraction of spectroscopic dropout sample showing strong Lya emission (WL ya > 50 A) as a function 
of UV luminosity for samples with 3.5 < z < 4.5 (blue circles) and 4.5 < z < 6.0 (red circles). Vertical lines correspond to ~ 90% 
completeness limits for B and V-drop samples. For the i'-drops, the completeness limits are ~ 0.7 mags brighter. Any incompleteness 
would serve to increase the Lya fractions further. 



the low end with Lya in the bluer filter). In principle, this 
could lead to Lya being preferentially recovered at higher 
redshifts. But we minimize these biases by simultaneously 
targeting the B, V, and i'-dropout population. For example, 
sources with very strong Lya emission at 5.5 < z < 5.7 (such 
that the line falls in the i^-band) may be scattered out of 
the z'-drop selection but would instead appear in V-drop se- 
lections. Thus, in this case, by conducting spectroscopy of 
V-drops, we can account for this diffusion of Lya sources. 
Similarly, our B-drop sample will contain the small num- 
ber of sources at 4.5 < z < 4.8 with very high strong Lya 
contaminating the V606-band filter (which would otherwise 
have little flux). While we don't target U-drops, Lya emis- 
sion from galaxies at 3.5 < z < 3.8 does not contaminate 
the £?435-band filter, so our B-drop sample should not have 
redshift-dependent biases. 

The redshift dependence of XL ya is affected not only 
by evolution in the internal properties of galaxies but also 
by the increase in the density of the IGM with redshift. 
At 2 ~ 6, the IGM provides a significantly greater optical 
depth to Lya photons than that at z ~ 4, resulting in a sec- 
ond order affect on the Lya fraction. In the absence of IGM 
density evolution, we would thus expect the redshift evolu- 
tion of the Lya fraction to be slightly greater than derived 
above. We can attempt to estimate the variation in XLya 
that is intrinsic to galaxy evolution (e.g. dust, ISM kinemat- 
ics) by subtracting the differential evolution expected from 
changes in IGM density. Deconvolving the effects of the IGM 
on Lya radiative transfer requires careful mode ling of the lo- 
cal density, velocity, an d ionisation field (e.g. ISantoj 12004 
iDiikstra et all l2007al lbl; IZheng et alj|2009h . We delay such 
a treatment to subsequent work s, and instea d we f ollow a 
very simple approach adopted in lOuchi et al.l (|2008l ) which 
yields a very rough estimate on the intrinsic redshift evo- 
lution of the Lya fraction. We compute the percentage of 
photons absorbed by the IGM assuming that the blue side 



of the Lya line is attenuated by exp[r a (z)], where T a (z) is 
the op tical depth for Lya photons as computed in iMeiksir] 
(2006). With this approach, we find that the IGM absorbs 
28, 42, and 49% of the Lya line flux at z ~ 4, 5, and 6. 
In a more sophisticated and realistic treatment, the den- 
sity and ionising background surrounding Lya emitters is 
likely to be greater than the mean, and infalling gas would 
also erase a fraction of the Lya line redward of rest-frame 
1216 A; the combination of these effects can cause the red- 
shift evolution in the transmission of Lya photons through a 
reionised IGM to be co nsiderably d i fferent than implied b y 
our model above (e.g.. ISantosll2004l ; IDiikstra et alj|2007bj ). 
We will model this effect in greater detail in the future. For 
the sake of clarity, here we define the intrinsic W L ya ,o,int as 
the rest-frame equivalent width that would have been ob- 
served if not for IGM attenuation, where the IGM absorp- 
tion is taken to follow the numbers derived above. Adopting 
a fixed intrinsic W L y a,o,int> we derive Lya fractions as above 
and find that the differential redshift evolution in XL ya in- 
creases to dxLy Q /dz ~ 0.12. In the next section, we attempt 
to understand the factors likely to be creating this redshift 
trend and the luminosity trend presented in the previous 
section. 



4.3 The factors governing the Lya escape fraction 

Earlier we demonstrated that prevalence of strong Lya emis- 
sion increases toward lower luminosities and higher redshifts. 
Here we discuss the factors that are likely governing the ob- 
served trends prior to exploring the use of higher redshift 
galaxies and their line emission as a probe of cosmic reion- 
isation. In §4.3.1, we examine whether trends in dust ob- 
scuration could potentially drive the observed Lya fraction 
relations. In §4.3.2 and §4.3.3, we discuss how the geomet- 
ric distribution and kinematics of the surrounding ISM may 
impact the Lya fraction trends we observe and discuss fu- 
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Figure 14. Average UV slopes of B-drops with and without 
strong Lya emission. Average UV slopes of B-drops (grouped 
in bins of Mpy) from the spectroscopic sample with strong Lya 
emission (WL ya > 50 A) are denoted by red triangles, while the 
average UV slopes rest of the B-drop sample is denoted by pur- 
ple circles. For reference, we overplot the average UV slope as a 
function of Mj/i/ for the entire p hotometric sample of B-drops in 
GOODS-N and GOODS-S from lStark et al.l j2009h as filled yel- 
low squares. The average UV slope from B-drops in the Hubble 
UDF is shown as open blue squares. Typical uncertainties in the 
UV slopes of individual galaxies are shown as error bars on the 
bottom of the plot; uncertainties in the mean take both this error 
and the number of sources in each bin into account. 



ture spectroscopic work in progress to provide further insight 
into the ISM of feeble high-z galaxies. Finally in §4.3.4 and 
§4.3.5, we close our discussion with a variety of factors possi- 
bly affecting Lya that are more difficult to directly constrain 
obser vat ionally. 



4-3. 1 Impact of dust extinction 

Previous observations have demonstrated that, among lu- 
minous LBGs, those objects showing Lya in emission tend 
to display b luer UV continuum slopes than those with Lya 



absorption JShaplev et al.ll200ll, 120031; IVanzella et al.ll2009l : 
IPentericci et al-lhoosl : iKornei et al.ll2009T l. While the pres- 



ence of dust may enhance Lya relative to the continuum 
|Neufeldlll99ll ; lHansen fc OhfeoOrj ; iFinkelstein et al.ll2008h . 
the results in Figure [T4l suggest that most often, the presence 
of significant quantities of dust generally leads to increased 
absorption of Lya photons relative to the continuum. 

Here, we examine whether similar trends are seen at 
lower luminosities. In order not to bias the UV colors, we 
determine the slopes using filters that are not contaminated 
with Lya emission or IGM absorption. Obtaining accurate 
UV slopes requires very accurate color measurements, thus 
we choose to focus on the B-drops, as the UV colors can be 
determined entirely from deep i' and z-band imaging with 
ACS. We translate t he UV colors into UV slopes using the 
relation presented in iBouwens et"afl (|2009bl ): /3=5.30(i 77 5 - 
2850 ) - 2.04. We have verified that this relation holds for a 
variety of star formation histories and age combinations that 
are appropriate for the z ~ 4 — 6 population. 

Given the relationship between UV color and j3, it is 



important to note that even a small photometric colour er- 
ror translates into substantial uncertainty in the derived UV 
slope. We thus estimate the luminosity limit at which the 
GOODS data become unreliable for estimating UV slopes 
by comparing the average UV slopes of B-drops in GOODS 
with UV slope measurements from higher S/N data in the 
UDF, using the photometric catalogs of ICoe et all (|2006f ) 
for our UDF sample. These results indicate that UV slopes 
measured from the GOODS dataset differ significantly from 
those determined from the high S/N UDF sample for B- 
drops with luminosities fainter than Muv — —20.5. There- 
fore, we cannot derive robust UV slopes for sources in our 
GOODS spectroscopic sample that are fainter than this 
limit. 

Concentrating on the brighter subset of objects, we ex- 
amine the UV slopes of our B-drop spectroscopic sample as 
a function of UV luminosity (Figure [T4")) . The data show that 
galaxies with strong Lya in emission (W hya , > 50 A) arc 
bluer than those systems without strong Lya emission, and 
are generally fit with f3 = — 2.0 across the ent ire luminosity 
range covered. Following iMeurer et all l| 19991 ). this value is 
consistent with very little dust obscuration (A1600 — 0.5), 
assuming a Calzetti extinction curve and normal stellar pop- 
ulations. The UV slopes of the overall population of B- 
drops are significantly redder than the Lya emitters, but 
when viewed as a function of luminosity, they grow steadily 
bluer toward low luminosities, ranging from f3 ~ —1.5 (at 
M UV ^ -21.5) to $ ~ -1.7 (at M UV ^ -20.5). The cor- 
relation between UV slope and emerging UV luminosity at 
these redshift s was first demonstrate d in a large photomet- 
ric sample in IBouwens et~ail J2009bl ): here we confirm this 
tr end with a spectroscop ic sample of dropouts. As argued 
in IBouwens et all l)2009bh . the trend most likely arises as a 
result of lower luminosity galaxies having less dust obscura- 
tion. Given that strong Lya emitters tend to arise in galaxies 
with little dust extinction, it is perhaps no surprise that we 
see a larger Lya fraction in low luminosity galaxies. Similar 
reasoning can also explain the observed redshift-dependence 
of the Lya fraction, as UV slopes a re found to grow stead ily 
bluer with redshift over 3 < 2 < 6 (|Bouwens et a l. 2009b). 

If we attribute the luminosity-dependence of the Lya 
fraction to the variation in dust obscuration, we can de- 
rive a relationship between mean UV slope and Lya frac- 
tion. From Figure [ill we find that d/3/dMuv — —0.23 over 
—21.5 < Afuv < —20.0. This is similar, but sligh tly steeper, 
than the value derived in IBouwens et all (|2009bV Over this 
same luminosity range, we derive dxLya/dMuv = 0.29. As- 
suming that dust extinction drives the evolution in the Lya 
fraction, we derive dxLya/dfi — —1.2. While this result was 
derived from the luminosity-dependence of XL ya , it should 
also apply to the redshift evolution if dust obscuration were 
to dominate the redshift-depend ent trends. Measureme nts of 
similarly luminous z ~ 6 LBGs (|Bouwens et al.ll2009bT ) sug- 
gest that the average UV slopes grow bluer by A/3 ~ —0.6 
from z ~ 4 to z ~ 6. If dust evolution dominates the redshift- 
dependence of the Lya fraction and the evolution in UV 
slopes, the above relationship suggests that x^ ya should in- 
crease by Ax Lya = 0.84. The fact that the actual redshift 
evolution is significantly less rapid (AxLya = 0.24 after 
roughly accounting for IGM absorption) can be explained 
by several factors. For example, it is likely that the lumi- 
nosity dependence of XLya is driven by more than just dust 
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obscuration (see the following sections for a discussion). Ad- 
ditionally as UV slopes grow bluer than /3 = —2.0, the varia- 
tion in the UV slope is possibl y driven by fa c tors other than 
dust o bscuration (e.g.. IStanwav et alJ feoOS: Bou wens et al.l 
2010a). Secondly, given that XLya cannot be larger than 1.0, 
the differential growth must slow down as XL ya increases; 
hence any extrapolation of the XL ya relations below our com- 
pleteness limits is highly uncertain. Finally, as emphasized 
in §4.2, our simple treatment of IGM absorption may in- 
correctly estimate the evolution in the Lya fraction that is 
intrinsic to galaxy processes. 

In summary, these results indicate that the variation 
in dust extinction with luminosity and redshift likely plays 
an important role in governing the observed evolution and 
luminosity-dependence in the Lya fractions. But, not sur- 
prisingly, it seems likely that additional factors play a role 
in governing XLya as well. We discuss these in more detail 
below. 

4-3.2 The hydrogen covering fraction 

After Lya photons escape the H II regions where they were 
created, they must diffuse through neutral gas and dust 
at larger radii. This could include gas that is participat- 
ing in outflows surrounding the galaxy in addition to gas 
that is being accreted onto the galaxies. Rest-UV spectra of 
high-redshift star-forming galaxies reveal low-ionisation ab- 
sorption hnes_tliat_areg_ej^aUv_bl^ ~ —200 
km s" 1 l|Shaplev et al.ll2003l ; Isteidel et alj|20ich relative to 
the centre of rest, indicating that outflowing neutral gas 
is nearly always present, while accreting gas appears less 
prevalent. The geometrical distribution and kinematics of 
this outflowing material plays a crucial role in governing the 
escape of Lya photons. This is clearly evidenced by the fact 
that Lya photons are typically observed to be redshifted 
by ~ 400 — 500 km s _1 with respect to the centre of rest 
(|Steidel et alJl201Ch . This result implies that the majority 
of Lya photons to escape through the absorbing material 
along the line of sight are those that achieve a sufficient (red- 
shifted) velocity such that they can travel through the inter- 
vening material without resonantly scattering. This is most 
easily accomplished in a model in which the Lya photons 
that are observed are those that have been "backscattered" 
off of the ouflowing material on the far side of the galaxy 
(e.g.. IShaplev et alJliooll : Isteidel et alJl20ld ). Hence, it is 
likely that the increased prevalence of Lya emitters among 
low luminosity galaxies tells us something about the distri- 
bution and/or kinematics of the ISM of feeble sources. In 
this section, we discuss the relationship of the ISM distribu- 
tion, dust obscuration, and Lya emission in luminous LBGs 
as found in previous work, and then consider whether a sim- 
ilar picture is likely to hold for low luminosity galaxies; in 
the following section, we focus on the kinematics of the ISM. 

Even for luminous sources, observational constraints on 
the distribution of the blueshifted neutral absorbing gas in 
the immediate vicinity of high redshift galaxies are lim- 
ited. The spectrosco pic study of z ~ 3 LBGs conducted in 
IShaplev et al.l (120031 ) revealed a strong correlation between 
the equivalent width of low-ionisation interstellar absorption 
lines (Wlis) and that of Lya, in the sense that the strongest 
Lya emitters tend to have the least absorption by the low- 
ionisation interstellar medium. As the absorption lines are 



highly saturated, IShaplev et all (|2003h note that the trend 
in Wlis is due either to variations in the velocity width or 
the covering fraction of absorbing gas, with the latter argued 
to be the dominant factor, such that sources showing strong 
Lya emission, on average, have the patchiest distribution of 
absorbing gas covering the continuum source (at least along 
the line of sight). Galaxies also s how a strong c orre lation 
between Wlis and E(B-V), which IShaplev et all l|2003h ar- 
gue implies a significant fraction of the dust which reddens 
the stellar continuum (and absorbs Lya photo ns) is located 
within the outflowing neutral gas. As noted in lSteidel et al.l 
(2010), from a theoretical standpoint, the exact effect of a 
non-uniform covering fraction on the observed Lya flux is 
not obvious. But regardless, these observations thus suggest 
a scenario in which strong Lya emission is generally coupled 
with low dust extinction and a low hydrogen covering frac- 
tion, a picture supporte d by recent observations o f strongly- 
lensed LBGs at z ~ 3 l|Quider et alj|2009l . I20T0I ) for which 
direct measures of the covering fractions of various ions are 
available. 

If similar trends are present in low luminosity galaxies, 
then the fact that Lya is much more common in feeble galax- 
ies may imply that these sources typically have lower cov- 
ering fractions of absorbing gas than more luminous LBGs. 
Addressing whether this is indeed the case requires deep 
spectroscopy of UV- faint systems, and thus is perhaps only 
feasible via studies of gravitationally-lensed galaxies. Some 
progress can be made with current field samples however. 
In particular, we can determine whether the coupling be- 
tween strong Lya emitters and weak ISM absorption is also 
present for feeble galaxies by creating composite spectra of 
Lya emitters binned by Muv- We will present the results 
of this analysis in a subsequent paper (Stark et al. 2010, 
in preparation). Some indication that the correlation is in 
place at low intrinsic luminosities is already apparent from 
the ^ -drop composite spectra presented in IVanzella et al.l 
(2009). As all the i'-drops in their sample are UV faint 
(< Muv >— —20), this spectrum provides insight into the 
properties of feeble sources. Owing to the faintness of the i! - 
drops, the composite is dominated by strong Lya emitters 
and shows very weak interstellar absorption lines. Higher 
S/N spectra are required to ensure that the absorption lines 
are saturated and to quantify the absorption line equivalent 
widths. 

Finally, while admittedly speculative, we note briefly 
that if a low hydrogen covering fraction is more common 
for low luminosity galaxies, it may also enable Lyman con- 
tinuum photons to more easily escape from feeble systems. 
Naively, this statement seems contradictory since the pres- 
ence of Lya photons stems from the absorption of ionising 
photons. However high-redshift galaxies are clearly not per- 
fect H II regions. In practice, a significant fraction of ionising 
photons are absorbed in the ionised regions surrounding the 
massive stars, leading to the production of Lya photons. 
Both Lya photons as well as any escaping ionising photons 
will then approach the surrounding neutral ISM, much of 
which is likely moving at great speeds with respect to the 
stars. Clearly, those systems with significant holes in their 
surrounding distribution of hydrogen will leak a larger frac- 
tion of ionising radiation. Indeed, the combination of strong 
Lya and significant ionisi ng photon escape fra ctions are seen 
in observations at z ~ 2 (Sh aplev et al.ll2006h . Intriguingly, 
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recent results reveal a possible trend toward greater Lyman 
continuum leakage in sources with low UV continuum lu- 
minosities (Steidel et al. 2010, in preparation). While much 
work is still required to verify the luminosity-dependence in 
the hydrogen covering fraction and the Lyman continuum es- 
cape fraction, these results strongly motivate detailed study 
of the physical properties of low luminosity galaxies. 



4.3.3 Kinematics of ISM 

As we discussed above, Lya is affected not only by the geom- 
etry and column density of the dust and hydrogen around 
it but also by the v elocity distribution of the surrounding 
neutral HI. Recently, ISteidel et all (|2010h has compared the 
kinematics of the absorbing gas (extracted from the proper- 
ties of rest-UV absorption lines) and the properties of Lya 
as a function of total baryon mass for a sample of z ~ 2 UV 
continuum selected galaxies. These results reveal that the 
low mass subset shows significantly stronger Lya emission, 
as we would expect from Figure [13] on the assumption that 
the feeble sources that we have studied spect r oscop ically are 
generally lower mass systems. ISteidel et al.l l|2010T ) demon- 
strate that the absorbing medium is significantly different 
for these two mass subsets, with the more massive galax- 
ies more often having a significant component of interstellar 
absorption at velocities close to the systemic redshift (and 
often extending to positive velocities). These high mass sys- 
tems not only have weaker Lya emission but the redshifted 
velocity of Lya is greater than in low mass galaxies. The 
nature of this excess absorption component at zero veloc- 
ity for massive gala xies is unclear; as discussed in detail in 
ISteidel et aD <|2010l > it could possibly be outflowing material 
which has stalled due perhaps to the larger gravitational 
potential, or alternatively infalling gas. 

Regardless, it is clear that there are significant mass- 
dependent variations in the kinematics of the absorbing gas, 
particularly at v = 0, such that lower mass galaxies transmit 
a larger fraction of their Lya photons to earth. These trends, 
together with those seen in dust obscuration and perhaps hy- 
drogen covering fraction, may play a large role in governing 
the luminosity-dependence of the Lya fraction seen in Figure 
1131 Future spectroscopic work of feeble galaxies which are 
magnified via gravitational lensing can help confirm these 
mass-dependent trends in the ISM kinematics. 



4-3.4 Addtional Factors Governing the Lya fraction 

We conclude by mentioning several additional factors which 
may contribute to the observed Lya fraction trends. Firstly, 
the Lya equiva lent width is highest i n the earliest stages o f 
star formation (|Charlot fc Fall 19931 ; iLeitherer et al.lll999h . 
when the contribution from massive, young stars is at its 
greatest. Studies of UV faint narrowband LAEs have gener- 
ally shown that these sources have younger ages and higher 
speci fic star formation rates than more luminous LBGs 
(e.g. iPirzkal et all 120071 ; lOno et al.ll2010h . However, a re- 
cent study of spectroscopically-confirmed LBGs has demon- 
strated that those sources with strong Lya emission tend 
to be older t han those without Lya emission at simi lar UV 
luminosities l|Shaplev et all 1200 ll ; iKornei et all 120091 ') . lead- 
ing the authors to suggest a physical picture whereby Lya 



emission escapes more easily once large-scale outflows have 
had time to reduce the dust covering fraction sufficiently. 
The sa mple presented in thi s paper allows us to extend the 
work of IKornei et al.l (|2009h to lower continuum luminosi- 
ties, which we plan to present a future study (Stark et al. 
2010, in preparation), taking care to account for the possi- 
ble c ontribution of nebular emiss ion lines to the mid-IR flux 
(e.g. ISchaerer fc de Barrosl [20091 ). 

Additionally, Lya equivalent width becomes larger at 
low metallicities owing t o the increase d ionising flux and 
harder UV spectra (e.g., [Schacrer 2003). Thus if low lumi- 
nosity galaxies have much lower metallicities than luminous 
systems, then we may expect to see an increased fraction 
of strong line emitters among UV faint systems. Only via 
future direct measurements of the metallicity of low lumi- 
nosity galaxies at these redshifts will we be able to predict 
the magnitude of this effect on the luminosity dependence 
of the Lya fraction. 

Finally, it is, in principle, conceivable that the trends 
could be driven by variation in the stellar IMF with lu- 
minosity an d redshift. As has been demonstrated in previ- 
ous studies l)Schaererll2003r ). top-heavy IMFs can boost the 
equivalent widths of Lya relative to that expected for nor- 
mal stellar populations. While several st udies have argued 
that the IMF may vary with redshift fe.g. lDavelkoOSr ). there 
is little direct evidence of such a variation. 



4.4 Comparison with narrowband Lya emitter 
studies 

Finally, we compare our results to those determined 
from studies of narrowband-selecte d Lya emitters (LAEs) 

~U [2006; 



in the same redshift range (e.g. IShimasaku et al.l _ 
iKashikawa et al.ll200rj : lOuchi et alll200gf ). We first contrast 
the two samples and examine whether the large fraction 
of strong Lya emission seen among our UV faint dropouts 
is consistent with the independently-determined luminosity 
functions of LAEs. Then we consider whether we can use 
the Lya trends discovered in our UV continuum samples to 
explain the lack of redshift evolution in the luminosity func- 
tions of the LAE population, which contrasts markedly with 
the strong evolution seen in the LBG samples. 

It is commonly asserted that narrowband LAE samples 
are fainter than UV continuum LBG samples. While this is 
true if one compares the LAE sample s to the well-studied 
spect roscopic z ~ 3 population (e.g. IShaplev et al. I l200ll . 
12003? ) . current photometric LBG samples extend to much 
fainter UV luminosities (e.g., Muv — —16 at z ~ 4 in the 
Ultra Deep Field) , and as we have discussed above, our spec- 
troscopic observations take advantage of these faint samples, 
extending to Muv — — 18 (Figure 2). Many LAEs identified 
in typical ground-based surveys are below the UV contin- 
uum magnitude limits and require stacking to determine the 
typical continuum lumin osity. In a recent analysis of a large 
sample of z = 3.1 LAEs. lOno et al.l (|2010h have shown that 
the stacked continuum magnitude is i' ~ 27, corresponding 
to Muv — —18.6, comparable to the UV continuum lumi- 
nosity of the faintest galaxies in our spectroscopic sample. 
Thus by probing down the LBG luminosity function in our 
DEIMOS survey, we are able to directly compare the Lya 
trends discovered in both populations. 

Consider the percentage of LBGs at a given redshift 
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that show strong Lya emission. Assuming that LBGs form 
the parent population of Lya emitters, we can compute 
the expected number density of strong Lya emitters as 
a function of Muv- Our spectroscopic data suggests that 
10% of B-drops with Muv — —21 show Lya emission with 
WLya,o > 50 A. Based on the number density of LBGs in 
this luminosity range, these data predict that the density of 
strong line emitters with Muv — —21 should be ~ 1.8x10 
Mpc -3 mag -1 . This value is very similar to the abundance 
of LAEs in this luminosity range from the UV luminosity 
functions of narrowband L ya emitters at z ~ 4 (1.2x10 
Mpc -3 ) (|Ouchi et al.ll2008T ). Clearly this calculation is very 
uncertain owing to the different selection and equivalent 
width limits, and is meant as a zeroet h order comparison . We 
note that of the 26 z = 3.7 LAEs in lOuchi et"aD l|2009h . 23 
have rest-frame equivalent widths greater than our threshold 
of ~ 50 A, so the LAE sample isn't probing significantly fur- 
ther down the equivalent width distribution then our LBG 
sample. At fainter continuum luminosities (Muv > —20), 
the agreement begins to break down, as our measured Lya 
fractions predict a larger ab undance of LAEs than inferred 
from the lOuchi et alj (|2008l ) U V luminosity functio ns. This 
is not surprising given that the lOuchi et alj (|2008T ) is only 
computed over —21.7 < Muv < —19.7, (where the data are 
sufficiently complete), requiring M uv and a to be fixed. 

Finally we consider whether the observed evolution of 
the Lya population over 3 < z < 6 is qualitatively consistent 
with the trends suggested by the luminosity and redshift- 
dependence of the Lya fraction. The rapid increase in the 
Lya fraction toward lower luminosities suggests that deep 
narrowband Lya samples will generally be dominated by 
UV faint galaxies. Since the redshift evolution of the num- 
ber density of UV faint LBGs is much less rapid than that 
of luminous LBGs, we expect the number density of Lya 
emitters (which should be weighted more toward UV faint 
sources) to decrease less rapidly than LBG samples. This 
trend is enhanced by the increase in the prevalence of Lya 
emission in LBGs over this redshift interval. Hence given the 
luminosity and redshift-dependent tren ds in the Lya frac - 
tion, it is not necessarily surprising that lOuchi et all (|2008T ) 
reveal that the observed Lya luminosity function does not 
evolve significantly with redshift over 3.1 < z < 5.7. 



5 IMPLICATIONS FOR REIONISATION 

In the previous section, we showed that strong Lya-emitters 
become more common between z ~ 3 and z ~ 6. We ar- 
gued that this trend is likely driven in part by a decrease 
in dust extinction and the covering fraction of hydrogen 
surrounding the H II regions where Lya photons are origi- 
nally produced. The first estimates of UV slopes at z > 6 
imply that the dust obscuration continues to decline to 
z ~ 7 — 8. Given these trends, the expected signal of reioni- 
sation (decreasing the prevalence of Lya emitters) should be 
readily apparent in deep spectroscopy of newly discovered 
zsso-band and Yips-band dropouts (e . g. [McLure et al.ll201o| ; 
Oesch et alll20ld;lBunker et alJl20l3 ; lBouwens et alj|2010al : 
Wilkins et al.ll2010al lbl). " 

Given the increase we find in the Lya-emitting frac- 
tion over 3 < z < 6 (Figure 1 13 p it is now interesting to 
examine whether there is any decline seen beyond z ~ 6 



such as might arise from an increasing neutral fraction in 
the IGM. Although our first set of z-drop observations were 
conducted in poor conditions, we can still provide prelim- 
inary constraints on xi ya if we adopt a sufficiently bright 
equivalent width and magnitude limit. The flux limits from 
the October 2009 run indicate that detection requires an 
equivalent width of at least 75 A (rest-frame) for sources in 
the luminosity bin -20.5 < Muv < -19.5 (§2.1). The Lya 
equivalent widths of sources fainter than this limit (all but 
two of the z-drops in the UDF) are not usefully constrained. 
Also, owing to the increased noise stemming from poor at- 
mospheric conditions it is not possible to detect Lya with 
typical equivalent widths in any of the sources studied at 
redshifts beyond z ~ 6.65. Over the redshift range in which 
we are sensitive to Lya emission, completeness simulations 
indicate that recovery rate of lines with Whyafi > 75 A is 
~ 50%. The photometric redshifts derived from the broad- 
band SEDs indicate that five of the eight z-drops are likely 
to lie at redshifts above z ~ 6.65. If this is true, these ob- 
jects would escape detection in our DEIMOS spectra. We 
take this possibility into account in our discussion of the 
Lya fractions below. 

Figure shows the overall evolution in the Lya frac- 
tion with redshift within this restricted luminosity range. 
The fraction of LBGs with Lya emission above our chosen 
threshold grows steadily over 4 < z < 6, reaching nearly 20% 
at z ~ 6. Based on our discussion in the previous section, 
we argue that this net change likely arises from a combi- 
nation of the redshift evolution in the dust and hydrogen 
covering fractions (which increases the Lya fraction toward 
higher redshift) and the IGM density (which decreases the 
Lya fraction toward higher redshift). 

If the line shown in Figure[8]is Lya, and if all z-drops in 
our restricted luminosity range ( — 20.5 < Muv < —19.5) lie 
at z < 6.65 (as is necessary for detection with DEIMOS 
in the October dataset), then our measurements would 
indicate that the fraction from the October run with 
WLya > 75 A (corrected for incompleteness) is marginally 
greater than that measured at z ~ 6, consistent with the 
slow increase over 4 < z < 6. However, given the very large 
errors in the Lya fraction, the z ~ 6.5 measurements are also 
fully consistent with a significant decline over 6.0 < z < 6.5. 
However we point out that the measured Lya fraction may 
actually be larger than we estimate here, for if five of the 
z-drops are indeed at z > 6.65, as suggested by their photo- 
metric redshifts, then the Lya fraction at z ~ 6.5 would be 
~ 2.6x larger than the datapoint in Figure [151 indicates. 
Clearly additional spectroscopy in improved atmospheric 
conditions is required for robust statements. 

Beyond z ~ 6.6, existing datasets imply that the Lya 
fraction may begin to decline. In contrast to the abundant 
samp les of robust LBGs n ow in place at z ~ 7 and beyond 
(e.g. iMcLure et all lioiol; iBunker et alj l20ld ; lOesch et ail 
2010; B ouwens et al.ll2010al ). searches for Lya emitters at 
z > 7 with narrowband filters have yet to reveal a large 
population of line emitters. Wh il e candidates have been 
identified (e.g. IStark et al.ll2007bl ; iHibon et ail 120091 '). con- 
firmation has proven challenging. Indeed, the highest red- 
shift Lya em itter with robus t spectroscopic verification lies 
at z = 6.96 l|lve et al.ll2006f ). Updated data from this sur- 
vey (Ota et al. 2008) suggest that the number density of 
Lya emitters at z = 7.0 is only 17% of that at z = 6.6. 
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Figure 15. Fraction of spectroscopic dropout sample with - 
20.5 < Mtjv < —19-5 showing strong Lya emission (Wi,y a > 
75 A) over 4 < z < 7. We estimate the Lya fraction at z ~ 6.5 
from a recent DEIMOS run targeting WFC3 z-drops. Owing to 
poor atmospheric conditions on this run, statistics are insufficient 
to robustly constrain the evolution in XL ya . However in contrast, 
the implied Lya fraction from the lack of Lya detection in a 
sample of 8 lensed z -dro ps studied wi t h nea r-IR spectroscopy in 
iRichard et~aH (2008) and lKneib et al". I j2004 l is shown by the up- 
per limit at z = 7.3. If these candidates are real and at high-z, 
then the implied Lya fraction at z ~ 7 is lower than measured 
at z ~ 6. Red triangles and associated errors correspond to mea- 
surements that will be possible if a large percentage of z-drops 
expected to be discovered with future observations (e.g. ULTRA- 
VISTA and WFC3 MOT) are followed up with spectroscopy. The 
quoted neutral fra ctions correspond to tho se implied from a map- 
ping derived from iMcQuinn et al.l (120071 ) and are only valid at 
z ~ 7. 



Attributing this decline in number density to a decrease in 
Lya transmission (ignoring the contribution from the de- 
clining galaxy number density), we would expect see strong 
differential evolution in the Lya fraction over 6.6 < z < 7.0. 
In principle, the DEIMOS observations should have easily 
detected such a decline in the initial sample of 17 z-drops 
observed in October 2009. Unfortunately, as discussed ear- 
lier, the poor conditions prohibit detecting Lya at z > 6.6 
for the faint sources considered in this work. Therefore we 
are unable to verify this result with our current DEIMOS 
dataset. 

However, the only detailed spectroscopic survey of z- 
drops currently in the literature l)Richard et al.l[200ot ) does 
offer tentative support for the downward trend in the Lya 
fraction. In their paper, seven gravitationally-lensed z-drops 
were observed with the NIRSPEC spectrograph on Keck 
l|McLean et al.l Il998l ) , and none showed evidence for Lya 
emission (with rest-frame equivalent width limits as low as 
~ 30 A). Likewise, a similar absence of Lya was seen in 
the z g 6.8 grav i tationally lensed g alaxy in Abell 2218 
jKneib et alj 120041 ; lEgami et alj 120051 ) . The typical lumi- 
nosities of these sources are considerably lower than the re- 
stricted range used in Figure 1151 hence one would expect 
them to have Lya fractions that are significantly greater 
than those observed at lower redshifts in this figure. In con- 
trast, we measure an upper limit to the z ~ 7 Lya fraction 
which is suggestive of a significant decline over 6 < z < 7. 



At the very large equivalent width thresholds adopted in 
this figure, it is unlikely that this decline could be pro- 
duced by incompleteness owing to the presence of OH emis- 
sion lines. While these NICMOS dropouts are generally not 
as robust as the z-dr ops being identified with WFC3 (see 
iBouwens et al.ll2009al ), the lack of Lya is tantalisin g in light 
of the narrowband L AE results l|Ota et al.l |200S| ) . Such a 
marked decline in the number density of Lya emitters should 
easily be detected in future follow-up of current WFC3 z- 
drop samples. 

We now explore the potential ramifications that these 
results may have for reionisation and discuss the possibilty of 
improving and extending the XL ya measurements in the fu- 
ture. We first discuss how one may use these measurements 
to quantify changes in the IGM. Various theoretical stud- 
ies have examined how an increase in the neutral fraction 
of the IGM impacts the transmission of Lya photons (e.g. , 
Santos! |2004 iFurlanetto et all l200rj ; IMcQuinn et~ai1 120071 ; 



Mesinger fc Furlanettoll2008 ). In this paper, we consider the 



models of IMcQuinn et alj (|2007l ). which make use of 200 
Mpc radiative transfer simulations to compute the effect of 
inhomogenous reionisation on the transmission of Lya pho- 
tons and the c lustering of Lya emitti ng galaxies. In Figure 4 
of their paper, IMcQuinn et alj (|2007l ) present the Lya lumi- 
nosity function at z = 6.6 for various HI fractions, ranging 
from xhi = 0.00 to xhi = 0.88. The greater the HI frac- 
tion, the more the amplitude of the Lya luminosity function 
decreases with respect to the fully-ionised case. The relative 
amplitude of the luminosity functions for various ionisation 
fractions is directly related to the decreased transmission of 
Lya photons associated with reionisation of the IGM. For 
example, the amplitude of the Lya luminosity function de- 
creases by a factor of ~ 120 when the HI fraction increases 
to 0.88. 

Our goal is to translate the measured Lya fractions into 
esti mates of the ioni s ation state of the IGM. We can convert 
the IMcQuinn et al.l (|2007l ) results into a mapping between 
Lya fraction and xhi as follows. First, we naively assume 
that in the absence of any change in the ionisation of the 
IGM, the Lya fraction would not evolve between z ~ 6 
and z — 7 (see below for an alternative approach). We then 
assume that the predicted decline in the amplitude of the 
Lya luminosity function (with respect to when the IGM 
is fully ionised) will result in an identical decrease in the 
Lya fraction. Hence if the HI fraction is 0.88 at z ~ 7, we 
should expect the Lya fraction to be 120 x lower than its 
measured value at z ~ 6, when the IGM is ionised. With 
these assumptions, the current measurements at z ~ 6.5 
indicate that the IGM is still largely ionised. 

However given that the Lya fraction evolves with red- 
shift in the absence of changes in the IGM ionisation state, 
our 'no-cvolution' assumption above is clearly simplistic. 
Indeed if the dust obscur ation continues to ev olve over 
z > 6 — 8 (as implied by IBouwens et ail l2010aT ). then in 
the absence of IGM ionisation state evolution, the Lya frac- 
tion is likely to be even larger at z ~ 6.5 — 7. If this is the 
case, then the mapping between Lya fraction and IGM ioni- 
sation state discussed above would be incorrect. To account 
for this, we extend the smooth evolution seen in the Lya 
fraction over 4 < z < 6 (when the IGM is highly ionised) to 
z ~ 6.5 and z ~ 7. We adopt this prediction as the baseline 
XLya value consistent with an ionised IGM. We associate 
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deviations below this value with evolution in the IGM ioni- 
sation state in a similar manne r as above. Fol l owin g this ap- 
pro ach, the lack of Lya in the lRichard et all l|2008h sample 
and iKneib et alj |2004) source suggests xhi > 0.3 at z ~ 7. 

With new ground-based surveys (e.g. UltraVISTA) and 
WFC3 programs set to reveal hundreds of z and Y-drops 
in the next several years, it should be possible to derive 
accurate Lya fractions over the redshift range 6 < z < 8. 
Indeed, we estimate that feasible campaigns with optical 
and near-IR multi-object spectrographs should easily be able 
to detect Lya equivalent widths as low as 30-50 A (rest- 
frame) for a sample of 200 z and Y-band dropouts with 
continuum magnitudes spanning 26 < m < 28, enabling 
much improved statistics on the fraction of LBGs with Lya 
emission and the possibility to measure Lya fractions as 
low as ~ 1%. With such a spectroscopic sample, it will be 
easy to identify whether the Lya fraction at 6.5 < z < 8 
shows a deficit with respect to measurements at z ~ 5 — 6. 
In fact, this dataset would enable the characterisation of 
Lya fractions more than 10 x lower than those expected at 
z ~ 7 in the absence of any change in the IGM ionisation 
state (Figure [15]). With these results, it should be trivial to 
identify the differential evolution in the Lya fraction that 
would be expe cted from reioni s ation ; in the context of the 
simulations of McQ umn et all (|2007f ). these results enable 
the IGM ionisation state to be tracked between fully ionised 
(xhi — 0.0) and almost completely neutral (xhi — 0.8). 



6 SUMMARY AND CONCLUSIONS 

In this paper, we have presented a new technique aimed 
at constraining reionisation. By measuring the luminosity- 
dependent fraction of continuum dropouts with strong Lya 
emission at z ~ 4, 5, 6, 7, and so on, we are sensi- 
tive to sudden changes in the transmission of Lya pho- 
tons that would be expected if the ionisation state of the 
IGM changes, complementary to past efforts to link evo- 
lution in the Lya lumino s ity function to reionisation (e.g. 
Malhotra fc Rhoads!l2004l ; iKashikawa et al.l [200^ ; live et all 
20061 : lOta et al.ll2008M ~The spectra not only provide infor- 
mation on Lya but also reveal constraints on the absorb- 
ing gas along the line of sight and provide a robust spec- 
troscopic sample for derivation of stellar populations (e.g. 
stellar masses, ages). Through this information, we seek to 
simultaneously provide insight into the evolution of the most 
important factors governing Lya transmission (dust obscu- 
ration, HI covering fraction and kinematics), and to thereby 
better isolate the impact of the IGM ionisation state on the 
evolving Lya fraction. 

To complete this Lya fraction test, we have con- 
ducted a large Keck/DEIMOS spectroscopic survey of 
dropouts in the GOODS-S and GOODS-N fields. The cur- 
rent Keck/DEIMOS sample includes 268 B-drops, 95 V- 
drops, 64 j'-drops, and 17 z-drops and uniquely extends to 
faint continuum magnitudes (Muv — —18). When combined 
with existing VLT/FORS2 spectroscopic observations of B, 
V, and i'-drops in the two GOODS fields, our sample con- 
tains 627 dropouts. We summarize our key findings below. 

1. We present tentative spectroscopic identification of 
a z-drop identified i n the GOODS ERS imaging CDF-S 
l|Wilkins et alj|2010a[ ). The galaxy shows a 7a emission fea- 



ture at z — 6.44, consistent with the photometric redshift 
expected from the broadband imaging. If confirmed, this 
would be the first WFC3 z-drop with spectroscopic identifi- 
cation in the literature. Future deep spectroscopic study of 
WFC3 dropouts should vastly increase the number of Lya 
emitters in this redshift regime. 

2. The prevalence of Lya emitters is greater among low 
luminosity galaxies. We find that the fraction of dropouts 
with Lya equivalent widths in excess of 50 A increases from 
~ 10% at Muv — -21.5 to > 40% for sources fainter than 
Muv — —19.5. This result is consistent with previous stud- 
ies comparing the UV and Lya lu minosity function fo r sam- 
ples of narrowband Lya emitters l|Ouchi et al J 12008). This 
finding demonstrates that a large fraction of galaxies at the 
faint-end of dropout samples are likely very similar to the 
high-EW and UV-faint LAEs found via narrowband filters. 

3. We find that the fraction of strong Lya emitters at 
fixed luminosity moderately increases toward earlier time in 
the redshift range 3 < z < 6. Binning our sample in two 
redshift ranges (3.0 < z < 4.5) and (4.5 < z < 6.0), we find 
that the fraction of Lya emitters with rest-frame equivalent 
width in excess of 50 A increases in each luminosity bin 
considered. 

4. We examined the possibility that these luminosity 
and redshift-dependent trends are driven by evolution in 
the dust obscuration of high-redshift galaxies. We find that 
the z ~ 4 sources with strong Lya emission have consider- 
ably bluer UV slopes, implying significantly less dust ex- 
tinction. This result is consistent with a range of obser- 
vational studies (e.g.. | Shaplev et ail 120031 ; iPentericci et al.l 
120081 : iKornei et alj [2009). In light of recent results revealing 
that dust obscuratio n decreases toward lower luminosities 
and hi gher redshifts (|Reddv fc Ste idcl 2009 ; iBouwens et al.l 
l2009bh , it appears that dust evolution plays a major role in 
governing the Lya fraction. 

5. We discuss the possibility that the covering fraction 
of hydrogen may be lower in low luminosity galaxies. In more 
luminous galaxies, strong Lya emission is generally cou- 
pled with low dus t extinction and low equivalent width ISM 
absorption lines (|Shaplev et al J [20031 ) . the latter of which 
appears to arise from a non-uniform covering fraction of 
neutral hydrogen jQuider et al.| [2009L Given the increased 
prevalance of strong Lya emitters and low dust obscuration 
at low luminosity, this result may suggest that galaxies with 
a non-uniform covering fraction may be much more com- 
mon among UV-faint systems. Strong Lya emitters with UV 
faint continua do appear to show weak low-ionisation ISM 
absorption (Stark et al. 2010, in preparation), but higher 
S/N and higher resolution spectra of gravitationally lensed 
galaxies are required to directly measure the covering frac- 
tion. We note that if the covering fraction of hydrogen is 
indeed lower for low luminosity galaxies, this would in turn 
imply that Lyman continuum photons may more easily es- 
cape from feeble galaxies, a trend that appears to be seen in 
z ~ 3 galaxies with deep Keck spectra (Steidel et al. 2010, 
in preparation). 

6. We have measured the fraction of z-drops with strong 
Lya emission in our spectroscopic sample, enabling us to 
constrain XL ya at z ~ 6.5. The estimated Lya fraction ap- 
pears consistent with the redshift trend seen over 3 < z < 6. 
However the current sample is far too small to rule out a de- 
cline to z ~ 6.5. Efforts to extend these studies to z ~ 7 have 
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thus far failed to locate Lya emission. The lack of strong 
Lya in the Keck spectra of 7 candidate LBG s in Richard 
et al. (2008) and the lensed z ~ 6.8 source in iKneib et al.l 
|2004 ) provides tantalising evidence that the Lya fraction 
declines by z ~ 7. T he availability of more robust WFC3 /IR 
z ~ 7 samples (e.g . lOesch et al.ll2010l ; iMcLure et aill2010l ; 
iBunker et al.l |2010| ) enables this possibility to be readily 
tested with the new generation of near-infrared multi-object 
spectrographs. 

7. Using recent simulations linking the evolution in the 
abundance of Lya e mitters to the IGM ionisation state 
(McQ uinn et al.l(2007l ). we examine the possiblity of placing 
constraints on reionisation with our spectroscopic sample. 
While we emphasize that the results are strongly model- 
dependent, if the evolution in the Lya fr action at z ~ 7 were 
to hold up, the iMcQuinn et alj l|2007h simulations would 
suggest that the neutral fraction of the IGM is xhi > 0.3 
at z ~ 7. More robust estimates of the z ~ 6.5 — 7 Lya 
fractions are needed to confirm these preliminary findings. 

The recent emergence of new photometric samples of 
dropouts at z > 7 with the WFC3 onboard HST provides a 
great deal of promise for extending this test in the future. 
Long exposures with upcoming near-infrared multi-object 
spectrographs will enable significant detections of strong 
Lya emission lines from sources as faint as J ~ 28, allowing 
the Lya fraction test to be extended to z ~ 8. 
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